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Abstract 

We present a study of semi-leptonic B — > Dlv decays in quenched lattice 
QCD through a calculation of the matrix element (L>|c7^6|S) on a 24 3 x 48 
lattice at (3 = 6.2, using an O (a)-improved fermion action. We perform 
the calculation for several values of the initial and final heavy-quark masses 
around the charm mass, and three values of the light-(anti)quark mass around 
the strange mass. Because the charm quark has a bare mass which is almost 
1/3 the inverse lattice spacing, we study the ensuing mass-dependent dis- 
cretization errors, and propose a procedure for subtracting at least some of 
them non-perturbatively. 
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We extract the form factors h + and h~ . After radiative corrections, we 
find that h + displays no dependence on the heavy-quark mass, enabling us 
to identify it with an Isgur-Wise function £. Interpolating the light-quark 
mass to that of the strange, we obtain an Isgur-Wise function relevant for 
B s — ► Ds^lv decays which has a slope —£' s = l^ + ^stat.) l^(syst.) at zero 
recoil. An extrapolation to a massless light quark enables us to obtain an 
Isgur-Wise function relevant for B — ► D^iP decays. This function has a slope 
~i'ud = 0.9^g(stat.) + 2 (syst.) at zero recoil. We observe a slight decrease in 
the magnitude of the central value of the slope as the mass of the light quark 
is reduced; given the errors, however, the significance of this observation is 
limited. 

We then use these functions, in conjunction with heavy-quark effective 
theory, to extract V& with no free parameters from the B — > D*£p decay 
rate measured by the ALEPH, ARGUS and CLEO collaborations. Using the 
CLEO data, for instance, we obtain \V ch \ = 0.037 + ] + l + t (tttsxttt;) ttt , 

' ' 1 c0 ' -1-2-1 Vl+/3 A i(l)y l+<5 Vm 2 ' 

where 5i/ m 2 is the power correction inversely proportional to the square of 
the charm quark mass, and (5 Al (l) is the relevant radiative correction at zero 
recoil. Here, the first set of errors is experimental, the second represents the 
statistical error and the third represents the systematic error in our evaluation 
of the Isgur-Wise function. We also use our Isgur-Wise functions and heavy- 
quark effective theory to calculate branching ratios for B^ — ► D^tv and 
S (s ) -»■ D^£P decays. 



PACS Numbers: 13.20.He, 12.15.Hh, 12.38.Gc, 12.39.Hg 

Key Words: Semi-Leptonic Decays of B Mesons, Determination of Kobayashi-Maskawa Matrix 
Elements (V c b), Lattice QCD Calculation, Heavy Quark Effective Theory. 



2 



I. INTRODUCTION 



Semi-leptonic decays of B mesons have been the focus of much activity in the last few 
years. Experimentally, their rather large branching ratios have allowed thorough studies of 
their properties. Theoretically, they have been a fertile ground for new ideas. Moreover, 
the interplay between these experimental studies and new theoretical ideas has led to a 
greater understanding of the flavour sector of the Standard Model and, in particular, to 
measurements of the less well-known Cabibbo-Kobayashi-Maskawa (CKM) matrix elements 
V cb and V uh 0. 

The main theoretical development in the study of hadrons containing a heavy quark, 
such as the b or c quarks, is undoubtedly the discovery of heavy-quark symmetry HH and 
the development of the heavy-quark effective theory (HQET) M , which describes the strong 
interactions of a heavy quark with gluons and light quarks at low energies. If one considers 
the masses of the b and c quarks to be much larger than the QCD scale, Aqcd ; one finds that 
the dynamics of the light quarks and gluons coupled to a b or a c quark become independent 
of this heavy quark's flavour and spin. In this limit, QCD exhibits a new 5'f/(4) S p inxflavour 
symmetry, known as heavy-quark symmetry, which acts on the multiplet (c |,c I, b |,6 [). 
This symmetry simplifies considerably the description of the decays of hadrons containing 
a heavy quark. For instance, the 20 form factors required to describe the semi-leptonic 
decays — > D^£u and B^ — > D^iu as well as the elastic form factors of B*f*) and 
mesons,0 can all be expressed in terms of two universal form factors, £ Uy d and £ s , known 
as Isgur-Wise functions [0], which parametrize the non-perturbative dynamics of the light 
degrees of freedom. £ Ut d describes the decays of mesons containing a heavy quark and a u 
or d antiquark, and £ s describes the decays of mesons containing a heavy quark and an s 
antiquark. Moreover, heavy-quark symmetry requires these Isgur-Wise functions to be 1 
when q 2 , the square of the four momentum transfer, is maximum [[|. 

In an earlier work ||, we obtained the Isgur-Wise functions £, u ^ from a lattice study of 
elastic D meson scattering. A similar approach, but with a different lattice action, was taken 
by Bernard et al. |J and led to very similar results. In the present paper, we extend our 
earlier work to include decays of the form P — > P'iv , where P^ is a heavy-light pseudoscalar 
meson composed of a heavy quark, with a mass around that of the charm quark, and 
a light antiquark, q. These processes are described by matrix elements of the vector current 
Q'^Q. These matrix elements can, in turn, be decomposed in terms of two form factors, 
h + (Lu;m Q ,m Q/ ) and h'(u]m Q ,m Q/ ), given by 

(p'(p')\QYQ\p(p)) . . Al > h+ , 



(v + v'yh m Q ,m Q i) + (v - v'yh (io;m Q ,m Q >) , (1) 



where v w = p (,) /M pW )U j = v-v' = (Mf, + Mp, - q 2 ) /2M P M P , and m QW is the mass of Q (,) . 

In the limit of exact heavy-quark symmetry, the two form factors become independent 
of the masses of the initial and final heavy quarks and 



1 The subscript s is used to distinguish mesons in which the light, spectator antiquark is s from 
those in which it is either u or d. 
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h (u;m Q ,m Q >) = 
h + (uj;m Q} m Q/ ) = , 



(2) 



where is an Isgur-Wise function of the type described above, whose exact functional 
form only depends on the quantum numbers of the light spectator antiquark. The only 
change we make to these quantum numbers in the present paper is to vary the light antiquark 
mass. For simplicity of notation, this dependence will be left implicit unless stated otherwise. 

For heavy quarks of finite mass, there are two sources of corrections to the simple results 
of Eq. (HD. The first is hard-gluon exchange between Q and Q' across the vector current 
vertex. The second results from the modifications of the vector current and meson states by 
higher-dimension operators in HQET. These latter corrections are proportional to inverse 
powers of the heavy-quark masses. Thus, we have 

h l (uj; m Q , m Q/ ) = (a 1 + (3 1 (uj; m Q , m Q >) + ^(uj; m Q , m Q >)) , (3) 

for i = +, — , where a + = 1, or = 0, (3 l represents the radiative corrections and 7*, the 
power corrections. It is important to note that these two corrections incorporate all of the 
mass dependence of the form factors h\ As defined in Eq. ([|), the Isgur-Wise function, 
£(u;), is renormalization-group invariant and normalized to one at zero recoil as required 
by heavy-quark symmetry Jl: 

m = 1 ■ (4) 

The radiative corrections can be evaluated analytically in QCD since they are perturba- 
tive. To quantify them, we use Neubert's short-distance expansion of heavy-quark currents 
0. He considers semi-leptonic B — > Div and B — > D*iv decays and computes radiative 
corrections to the corresponding heavy-quark matrix elements to order a s as a function of 
m c and mj. His calculation improves the previous leading logarithmic evaluation of these 
corrections || in two ways: firstly, he includes next-to-leading logarithms in running the 
O ((m c /mb)°) heavy-quark operators from nib down to scales at which HQET can be safely 
used, and secondly, he obtains, to order a s , the full dependence of the heavy-quark current 
on the mass ratio z = m c /mb. The sum of these new contributions is as large as the lead- 
ing logarithmic term. Corrections to Neubert's computationi are of order a 2 s (z In z) n with 
n = 0,1,2 and should be smaller than 1%. The fact that Neubert's result accounts for the 
full order a s dependence of the heavy-quark current on the mass ratio z is important for us, 
because our range of heavy-quark masses is quite small (see Table |): z ranges from 0.6 to 
1. 

The power corrections are proportional to powers of egw = A/(2mg(;)) where A is the 
energy carried by the light degrees of freedom in the mesons. A will of course depend on what 
these light degrees of freedom are. In what follows, we will use A = A x = 500MeV when 
working with light degrees of freedom with spin 1/2 and isospin 1/2. Because €q(o ~ 1/6 for 
the heavy quarks we are considering, we would naively expect power corrections in h + (u)) 
and h~(oj) to be of order 15 to 30%. These corrections are difficult to quantify because they 



2 Neubert runs the O (m c /rrib) contribution at one loop. 
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involve the light degrees of freedom and are therefore non-perturbative. Luke's theorem ||10|| . 
however, guarantees that there are no O (eg) corrections to h + (u) at zero recoil and one 
may expect that power corrections to h + remain small away from zero recoil. This is not 
expected to be true for h~ which is not protected by Luke's theorem. 

For degenerate transitions where Q = Q', conservation of the vector current Q^Q 
provides further constraints on the radiative and power corrections: 

(3 + (l;m Q ,m Q ) = 
7 + (l; m Q ,m Q ) = 

l3-(uj;m Q ,m Q ) = (5) 
l~(u;m Q ,m Q ) = , 

where the last two equations hold for all uj. 

Our results come from a quenched simulation on a 24 3 x48 lattice at (3 = 6.2 on a sample 
of 60 gauge field configurations Jll]. The lattice has an inverse lattice spacing of around 



2.7 GeV \T2\. We do not suffer much here from errors associated with uncertainties in the 
determination of the lattice spacing since our main results are dimensionless and depend at 
most logarithmically on the scale. Our light quarks have masses which bracket the strange 
quark mass. Because our heavy quarks have masses in the region of the charm-quark mass 
which are large in lattice units (up to a half or more), we must contend with discretization 
errors proportional to powers of cmtiq, where tuq is the mass of the heavy quark. In order 
to reduce these discretization errors, we use the (9(a)-improved fermion action originally 



proposed by Sheikholeslami and Wohlert [13 with which discretization errors in operator 



matrix elements and hence in our form factors are reduced from O (arriQ) to O (a s amo) |L4 

The remainder of the paper is organized as follows. In Section [TT] we present the details 
of our simulation, as well as our strategy for obtaining the form factors h + and h~ from the 
calculated three-point functions. In Section [TT1] we discuss discretization errors and describe 
a procedure which enables us to subtract some of these errors non-perturbatively. In Section 
|TV] we present our results for the form factors h + and h~ for three values of the light antiquark 
mass and all available initial and final heavy-quark combinations. We also extrapolate h + 
in the light-antiquark mass to the chiral limit, and interpolate it to the strange quark mass. 
In Section [V] we study the dependence of h + and h~ on heavy-quark mass and attempt to 
extract the leading power corrections. We find that h + displays no measurable dependence 
on heavy-quark mass which enables us to conclude that this form factor is an Isgur-Wise 
function once radiative corrections are subtracted. In Section [VI], we study the dependence 
of h + on the light-quark mass and extract the Isgur-Wise functions £, u ^ and £ s . We find that 
the slopes of these functions at u — 1 are 

C(l) = -[0.9l^(stat.)^(syst.)] (6) 

and 

^(l) = -[l.2^(stat.)+?(syst.)] . (7) 

We thus observe a slight decrease in the magnitude of the slope light-antiquark mass; given 
the errors, however, the significance of this observation is limited. We compare our results 
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for these Isgur-Wise functions to other theoretical as well as experimental determinations. 
We find excellent agreement with experiment. In Section [VI]| we use our Isgur-Wise function 
£ M d to extract the CKM matrix element from different experimental measurements of 
the differential decay rate for B — > D*£u decays. Our results for \V c b\ are summarized 
in Table |X-V11| and are compared to other determinations of this matrix element. Our 



procedure for extracting \V c b\ differs from that proposed by Neubert |L5] in that we fix 
the uj dependence of the differential decay rate with our calculation instead of obtaining it 
from experiment. This enables us not only to extract |V^&| with no free parameters, but 
also to check the validity of non-perturbative QCD against experiment. We find that the uj 



dependence predicted by our calculation agrees very well with the results of the ALEPH [16 



and CLEO |L7| collaborations. In Section [VII]] , we use and £ s to compute the branching 



ratios for — > D^iu and B^ — > Dfyiu decays, and our results are summarized in 

d- 

3.2 + ^(lat.) ± l.O(hqs) (8) 



Table |XVIII| . We also compute ratios of semi-leptonic widths and find 

F(B -> D*£u) 



T(B — > Div 
and 

T(B S -> D\lv 



T(B s ^D s £u] 



3.3i?(lat.) ±1.0(hqs) , (9) 



where the first set of errors was obtained by adding our lattice statistical and systematic 
errors in quadrature and the second set of errors, denoted by "hqs" , quantifies the uncertainty 
due to neglected power and radiative corrections. We confront our predictions for these 
branching ratios and ratios of widths with experimental measurements where available and 



find that they compare quite favourably. Finally, in Section IX we present our conclusions 



II. DETAILS OF THE CALCULATION 

A. Lattice Action and Operators 

Since we are studying the decays of quarks whose masses are large in lattice units, we 
must control discretization errors. In order to reduce these errors, we use an O (a)-improved 
fermion action originally proposed by Sheikholeslami and Wohlert |13||, given by 



S S F W = Sj-i^Yl q{x)F iat {x)a iat q{x), (10) 



where S]¥ is the Wilson action: 
S f = E j 2(z)?O0 -K£[g(2)(l - 7 M )C/ M O0g(z + £) + q(x + /*)(! + l^U^qi 



x 



The leading discretization errors in matrix elements for heavy-quark decays obtained from 
numerical simulations with the fermion action Eq. ([H]) are reduced from O (arriQ) to 
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O (a s amq) and O (a 2 mgj, provided one also uses "improved" operators obtained by "ro- 
tating" the field of the heavy quark, Q: 

Q( x )^(l-l T D)Q(x) . (12) 

Thus, to obtain an O (a)-improved evaluation of the matrix element of Eq. ([!]), we use a 
"rotated" vector current 

Vt = Q\x)t^Q{x) , (13) 

where 

r^ = (1 + ^7" 5) 7^(1-^7-5) (14) 
and where the subscript / indicates that Vj is an improved lattice current. 



B. Extended Interpolating Operators 

In order to isolate the ground state in correlation functions effectively, it is useful to 
use extended (or "smeared") interpolating operators for the mesons. In this study we use 
gauge-invariant Jacobi smearing on the heavy-quark field (described in detail in Ref. |P-8||), 
in which the smeared field, Q s (x,t), is defined by 

g 5 ( x ,t) = ^K(x,x')g(x',t), (15) 



where 



and 



N 

K(x,x') = i^A n (x,x') (16) 

n=0 



A(x, x') = ^{5 xV _^(x - i, t) + 5 xV+ ^(x, *)}. (17) 
i=i 



Following the discussion in Ref. [fig], we choose ks = 0.25 and use the parameter iV to 



control the smearing radius, defined by 



r 2 = 



Sx|x| 2 |K(x,0)| : 
E x \K(x,0)\ 2 



(18) 



We use iV = 75, giving r « 5.2. 

In terms of the operator Q s of Eq. ([T5|) , the spatially extended source Jp we use to create 
pseudoscalar mesons composed of a heavy quark Q and a light antiquark q is given by 

J P (x) = q(x)(l + l - T D) 7 5 (1 - \r D)Q s {x) . (19) 
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C. Three- Point Functions and Lattice Form Factors 



The computation of the matrix element (P'(p') |Q'7 M Q|P(p)) proceeds along lines similar 
to earlier calculations of the electromagnetic form factor of the pion and to determinations 
of the form factors corresponding to semi-leptonic decays of the D meson into light mesons. 
(For recent reviews of lattice computations of weak matrix elements and references to the 
original literature see, for example, the reviews in Ref. ||19|| ). Thus, we calculate the three- 
point correlator, 

Cftfjp'.qW = £ e^e*'*(J P ,{t fi y)Vf{t,TL) 4(0,0)) , (20) 

where Jp is the spatially-extended interpolating field for P defined in Eq. (0), Vj 1 is the 
C(a)-improved vector current of Eq. fllcf) and p = q + p'. To evaluate these correlators, we 
use the standard source method reviewed in Ref. |2Tfl . 

Provided the three points in the correlator of Eq. (p0| ) are sufficiently separated in time, 
the ground state contribution dominates and 

<tf(f;p',qW -+ M r2lT (p/2) e- E *Wf-* (P'(p')\Vf(0)\P(p)) , (21) 

where Ep (Ep/) is the energy of the initial (final) meson and Zp(p 2 ) is the matrix element 
(0| Jp(0)|P(p)). To cancel the above time-dependence, we normalize the three-point function 
by two two-point functions and consider the ratio 

R(t - p -' l)Q ^' = c^, P ) Q c 2 (t,- t , P % / - (22) 



where 



C a (f,p) Q = Ee- !p ' x (Jp(«,x)4(0)> (23) 



and 



C 2 (t, p) Q e ' EpV2 cosh (^( T / 2 " *)) • ( 24 ) 

Here, T is the temporal extent of the lattice. (For t <C T/2, exp(— E P T/2)cosh(Ep(T/2 — t)) 
— > (l/2)exp(— Ept)). Thus, in terms of the form factors defined in Eq. (j^) 



R»(t;p',q) Q ^ Q i — + - ^MpM' P 

t,t f -t->oo Zp(p^)Zpt(p' z J v 



x ((v + vYh^ t (uj;mQ,mQ,) + (v - u') M /i lat .(^; Wq, m Q 0) , (25) 

where hf at are related to the continuum form factors, h ± , by a multiplicative renormalization, 
up to discretization errors, as discussed in Section |TJ. 
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To obtain the desired form factors, we fit the ratio W" of Eq. fl22|) to the asymptotic form 
of Eq. (|25D by minimizing, with respect to the parameters h^ t and /tj~ t , a x 2 function which 
takes into account correlations between the different times (labelled by t), but not between 
the different equations (labelled by /i). We neglect correlations between equations, because 
spatial and temporal components of Eq. (^) may be affected differently by discretization 
errors, as we discuss at the end of Section [III Q . The x 2 value that we quote indicates not 
only whether our ratios are asymptotic, but also whether the decomposition of in 
terms of h^ t (uj) and /ij~ t is good. In fitting the ratio we fix the wavefunction factors 
Z P (i) , the energies, E P (,) , and masses, M P ^ , of the mesons to the values obtained from a fit 
of the relevant two-point functions to the asymptotic form of Eq. (|24|), taking into account 
correlations in time. 

We first obtain hf^ (u) from the time component of Eq. (|25| ) alone, assuming that the 
contribution proportional to h^iu) can be neglected. This approximation is exact, up 
to discretization errors, for degenerate transitions, i.e. transitions in which the initial and 
final heavy-mesons are the same, and true up to radiative and power corrections for non- 
degenerate transitions, i.e. transitions between mesons which contain the same light anti- 
quark, but different heavy quarks (see Eq. ([5])). For these non-degenerate transitions we can 
get a posteriori some idea of the size of the contribution of h^ t (to) to the time component 
of Eq. fl25"D . Holding h^ t (uj) fixed to its time-component value, we use all non- vanishing 
components of Eq. ( |25"1) to obtain h{ at (u). We find (see Section |IV A|) that h^ t .(^)'s con- 
tribution to the time-component of Eq. (|25|) is less than about 1%, thereby justifying the 
approximation we make in obtaining h + (uj)i at 



D. Lattice Parameters and Details of the Analysis 



We compute the three-point function of Eq. ( p20|) for four values of both the initial and 
final heavy-quark hopping parameters, kq and Kq> taken from 0.121, 0.125, 0.129, 0.133 
(see Table 0]); three values of the light antiquark hopping parameter, n q (0.14144, 0.14226, 
0.14262); two values of the initial meson momentum ((0,0,0) and (1,0,0) in lattice units); 
and ten values of the momentum carried by the vector current (qa(12/7r) = (0, 0, 0), (1,0,0), 
(0,1,0), (0,0,1), (1,1,0), (1,0,1), (0,1,1), (1,-1,0), (-1,0,1), (0,1,-1)). To improve statistics, 
we average the ratios of Eq. (^2|) over all equivalent momenta. Moreover, data with initial 
or final momenta greater than 7r/12a are excluded because they have larger systematic and 
statistical uncertainties. Finally, we choose tf, the time at which the final meson is destroyed 
(see Eq. (p0|)), to be half-way across the lattice (i.e. tf = 24) and symmetrize the three- 
point functions about that point using Euclidean time reversal, also to reduce the statistical 
errors. 

We observe a plateau in the ratio R^it) of Eq. (|22D around t = 12, typically extending 
over 5 time slices. Therefore, we fit the ratio R^(t) over the range t = 11, 12, 13 to the form 
given in Eq. (|25"D for all momentum and heavy-quark mass combinations. For the purpose 
of illustration we plot, in Fig. [I], the ratio R°(t) vs. t for the case where the initial meson 
has momentum (7r/12a, 0, 0) and the final meson, momentum (0, 0, 0). We fit the two-point 
functions to the asymptotic form of Eq. fl2"4] ) in the range t — 11 to t — 22. The results of 
these later fits are given in Table |I[ 

Statistical errors are obtained from a bootstrap procedure [2~I]. This involves the creation 
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of 200 bootstrap samples from the original set of 60 configurations by randomly selecting 60 
configurations per sample (with replacement). Statistical errors are then obtained from the 
central 68% of the corresponding bootstrap distributions as detailed in Ref. [Illl . 



Use of the HQET implies a choice of the expansion parameter, mq, and this requires 
some care |22| , |2"3|j . We define trq as follows: 

m Q = ^(3M* + M*)-A x , (26) 

where Mp and M v are the relevant, chirally extrapolated pseudoscalar meson and vector 
meson masses in lattice units (see Table |). Since these masses correspond to heavy-light 
mesons whose antiquark is massless, the light degrees of freedom carry an energy A x = 
0.50 GeV as discussed after Eq. (0). 

In Table [TIT] and Table [TV| , we tabulate the results that we obtain for the radiative 
corrections, f3 + {u); ttiq, m,Qi) and (3~{u\ rriQ, mg/), of Eq. (|3|) for various combinations of the 
heavy-quark masses and for a few values of uj. As mentioned in Section |, we determine 
these corrections with the help of Neubert's work 0. Since our results for the form factors 
are obtained in the quenched approximation, we set the number of quark flavours to zero 
and assume no particle thresholds in Neubert's expressions!. 



III. Z v , DISCRETIZATION ERRORS AND HOW TO SUBTRACT THEM 

NON-PERTURBATIVELY 

Throughout this study of semi-leptonic weak decays of heavy mesons, we use an O (a)- 
improved fermion action and take for the lattice vector current, the "improved" operator 
Vf 1 of Eq. (|13"D, as discussed in Section [II A| . In concrete terms, this means that we expect 



mass-dependent discretization errors to be of O [a s amQ) ~ 5% and O ((amQ) 2 ) ~ 10% at 
the charm mass I instead of O (amq) ~ 40% and O ((amg) 2 ) ~ 10% as they would be 
without O (a)-improvement. Thus, despite the improvement expected, discretization errors 
in our calculation could be significant. 

Discretization errors in the lattice evaluation of the matrix element of Eq. ([[]) can be 
parametrized as follows: 

(py)IQ'r>Q|p(p)) lat , 

Z v (a s (a)) . = 

(v + v'f (l + d + (uj)) h + (uj) + (v- v'Y (l + dT{u)) h~(oj) + O (a 2 ) , (27) 



3 There is, in fact, no rigorous way of running quenched lattice QCD results since the lattice cutoff 
a -1 is adjusted to incorporate in part the effects of quenching. 

4 For this estimate, we use the boosted value of the coupling constant g 2 oost = (8K cr it) g 2 ^ 1.66 
and the improved bare mass defined before Eq. (33) with kq = 0.129. 
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where Zy{a s (a)) is the usual renormalization constant which relates the lattice vector cur- 
rent to the continuum one @. Because this constant describes physics that takes place above 
and around the lattice cutoff, it is perturbative and independent of the initial and final 
states. 

In Eq. (|27D d + and d~ are the Euclidean-invariant discretization errors to all orders 
in a. At O (a 2 ) the hypercubic group allows for additional errors which depend on the 
Lorentz index of the vector current. The discretization errors are non-perturbative and 
depend on the initial and final states, because they correspond to matrix elements of higher- 
dimension operators which are artefacts of lattice regularization. In addition, they depend 
on the procedure used to cancel all the factors which relate the three-point function to the 
matrix element (see Eq. (pT|) ) . We adopt the expedient of assuming that we can absorb the 
Euclidean-invariant discretization errors into an effective renormalization constant Zy. 

In the remainder of this section, we will attempt to quantify the discretization errors in 
our calculation more precisely and describe a procedure which enables us to subtract them, 
at least partially. 



Determination of Zy* 



To study discretization errors, we define an effective renormalization constant, Z v - , for 
vector currents composed of degenerate quark fields (i.e. of the form q^^q) by 

7 eff _ 1 C 2 (tf,p) , , 

for tf = T/2, where T is the temporal extent of the lattice, C3 1 and tf are defined in 
Eq. (|20|) and C2 in Eq. (|23|) . In the absence of discretization errors, Eq. ( p8|) yields a very 
accurate non-perturbative determination of the renormalization constant Zy. To see that 
the ratio of Eq. (|28|) is in effect Zy, one must use the fact that the forward matrix element 
of the temporal component of the vector current is the charge, up to a trivial normalization 
factor. The factor of 1/2 comes from our boundary conditions (see Eq. fl24|)). Unless stated 
otherwise, we will take p = 0. In the presence of the discretization errors described in 
Eq. (pTD, however, the ratio of Eq. (^) becomes 

Z e v ff = Z v (l-d + (l) + 0(a 2 )) . (29) 



We start the discussion with a review of the determinations of Zy for currents composed 
of degenerate light-quark fields, between pseudoscalar states composed of degenerate, light 
quarks and antiquarks where we expect Zy* to be close to Zy. Using 10 gluon configurations 
from our simulation at (3 = 6.2, we find 



5 It is important to note that similar discretization errors are present for all definitions of the 
current, even the conserved current away from the forward direction. 
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Z$ f = 0.8314(4) at « = 0.14144 
Zy ff = 0.8245(4) at « = 0.14226 

= 0.8214(6) at « = 0.14262 (30) 

These results confirm that discretization errors are small for light quarks (less than about 
2%), and we take 

Z v = 0.82(1) (31) 

as our best estimate for Zy. This value is also consistent with the expectations from one-loop 
perturbation theory [p5[| : 

Z v = 1 - O.Wg 2 + 0(g 4 ) ~ 0.83 at /3 = 6.2 (32) 

when evaluated using the boosted value of the coupling constant, obtained from the mean 



field resummation of tadpole diagrams ]26[ . 

We now turn to the evaluation of Zy using Eq. for degenerate heavy-quark currents 
between pseudoscalar mesons consisting of a heavy quark Q and a light antiquark q. The 
results and, in particular, the difference from the value in Eq. (|3"1~D, give us a measure of the 
size of the discretization errors, which are of O (« s flmg) and O (a 2 mq^ here. In Table [V], 

we present the results for Zy, obtained from the simulation at /3 — 6.2 for four values of 
the heavy-quark mass, and with the light-quark mass corresponding to K q = 0.14144, and 
from a simulation at (3 = 6.0 for three values of the heavy-quark mass and with the light 
quark hopping parameter equal to 0.144 i. Also tabulated are estimates of the improved, 
bare mass of the heavy quark, rag, defined by amq = amn(l — (l/2)amg), where amq = 

(1/2)(1/« Q -1/K 0rft ). 

In Fig. ||we plot the results for Zy* as a function of rriqa for the two values of (3. Fitting 
this behavior to a quadratic function of rriqa, 

Zy ff (K Q ) =A + Bm Q a + C{m Q af (33) 

we find A = 0.814(2) (A = 0.791(4)), B = 0.342(12) (B = 0.397(18)) and C = -0.072(18) 
(C = -0.120(20)) at P = 6.2 ((3 = 6.0). These fits are excellent. It is interesting to note that 
the results extrapolate to approximately 0.81 (0.79) in the chiral limit and are thus in good 
agreement with the values determined using light quarks as can be seen in Fig. |2| where we 
have also plotted the light-quark values for Zy* given in Eq. (|30D . This fact together with 
the observation that the size of the mass-dependent effects for a given arriq is very similar 
at the two values of (3 gives us confidence that the mass dependence we observe is indeed 
due to discretization errors. 



Further results from the simulation at (3 = 6.2 are presented in Table and in Fig. y. 
For kq = 0.129 and 0.121 we have evaluated Zy 3 F (kq) at three values of the mass of the 



6 The simulation at j3 = 6.0 was performed with 36 quenched gauge field configurations on a 
16 3 x 48 lattice using the O (a)-improved SW action of Eq. (|To|). For details of the simulation, 
please see Ref. [27 
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light quark. The results can be seen to be practically independent of the mass of the light 
quark. We have also evaluated Zy^ using Eq. (H) with p = (vr/12, 0, 0) and n q = 0.141440. 
The difference between the results obtained with p = (7r/12a,0,0) and with p = is less 
than 1%. Finally, we have determined Zy^Kq) using 

for p = (7r/12a, 0, 0) and tf = T/2 and where Ep(p 2 ) is the energy of the meson with 
momentum p. Now it is no longer the charge operator which appears in C3, and the 
statistical errors increase significantly (see Fig. [|). The values of Zy^ given by Eq. fl3"l"D are 
consistent with those obtained with \i — to within 1.5 standard deviations. 



B. Implications of the Results for Zy 

The results for Z v (kq) with p = presented above differ from the value of Zy given 
in Eq. (|31~D by about 10-20% for the range of quark masses used in our simulations (for 
kq = 0.129, which corresponds approximately to the charm quark for both values of /3, the 
difference is about 12% ). This difference is a good indication of the size of mass-dependent 
discretization errors in our calculation; it is consistent with our expectation that they should 
be of O (Q s flmg) and O (a 2 m;Q. 

Our results for Zy also enable us to quantify the dependence of discretization errors on 
momentum as well as on the Lorentz component of the current used to obtain them. As noted 
in the previous subsection, the difference between the results obtained with p = (7r/12a, 0, 0) 
and with p = is less than 1%. This is a clear indication that as long as we limit ourselves 
to momenta p such that |p| < n/12a, discretization errors proportional to ap are small. 
As for the dependence of Zy on the Lorentz index of the current, the situation is less 
clear. The ratio Z v ff (0.121; n = l)/^ (0.121; fi = 0) for p = (7r/12a,0,0) indicates 
that this dependence could be as large as 11%. However, given that the statistical errors on 
Zy (0.121; fj, — 1) are quite large, much of this dependence could be a statistical fluctuation. 



C. Non-Perturbative Subtraction of arriQ-Errors 

Having isolated and quantified the different sources of discretization errors, we now 
investigate the possibility of subtracting these errors. It is important to remember that 
these discretization errors are given by matrix elements of higher dimension operators: they 



7 The statistical errors in Zy* (kq) are tiny, due to a cancellation in the ratio ( p8| ) of the fluc- 
tuations in the numerator and denominator. In order to get such a dramatic cancellation of the 
fluctuations it is necesary to have precisely the same momentum in the numerator and denomina- 
tor. If, for example, we take p = (ir/12a, 0,0) in C3 but average over all 6 equivalent momenta in 
C2, ((±7r/12a, 0, 0), (0, ±7r/12a, 0), (0, 0, ±7r/12a)), then the statistical error in the ratio increases 
enormously. 
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are non-perturbative and will depend on the initial and final states between which the current 
Vf is sandwiched. This means that in any attempt to subtract them, one must evaluate 
the relevant corrections with states as similar as possible to the ones which appear in the 
matrix element of interest. With this in mind, we have devised the following subtraction 
procedure. 



Firstly, as mentioned in Section III, we assume that the mass-dependent discretization 



errors can be absorbed into an overall effective normalization: 



<p'( P ')io'f^iP(p)) la , = ^r/'r/ , (35) 



(P'(p')\Q'rQ\P(v)) 

Zy ff (aMp,aM P ,]fi) 

where T M is defined in Eq. (|T4]). 

Secondly we find a normalization condition, i.e. a kinematical point at which we 
know the physical value of the matrix element. For the case of degenerate transitions, 
this normalization condition is simple; electromagnetic charge conservation requires that 
h + (l;mQ,mQ) = 1. For the case of non-degenerate transitions, the normalization condition 
is slightly more complicated. HQET requires, as we saw earlier, that h + (l;mQ,mQi) — 1 + 
nig, nig/) + 7 + (l; m<g, m<g/). The radiative corrections, toq, rriQi), we know from 

perturbation theory. The power corrections, 7 + (l;mQ,mQ/), are non-perturbative and are 
yet to be determined in a model-independent and reliable way. We are, however, helped here 
by Luke's theorem which guarantees that h + (l; ttlq, irqi) is free of corrections proportional to 
a single power of the inverse heavy-quark masses. Thus, 7 + (l; tuq, mqi) ~ e %,Q> + @ ( e Q,Q') 
and is small. In fact, as we shall see shortly, the exact size of 7 + (l; rriQ, rriQi) is not important 
for determining the Isgur-Wise function. Thus, we will take our normalization condition to 
be 

h + (l;m Ql m Q ,) = 1 + i3 + (l;m Q ,m QI ) (36) 

for both degenerate and non-degenerate transitions. 

This condition determines Zy. With Zy* defined by Eq. fl3"5|) we find 



Kd l \ m Q, m Q') V ' 

where hf &t (l;mQ,mQ>) is the zero-recoil form factor obtained from our lattice calculation 
and the O (a 2 ) stands for discretization errors which are not Euclidean invariant. Because, 
as we mentioned earlier, discretization errors made in the evaluation of a matrix element 
depend not only on the initial and final states considered, but also on the procedure used to 
obtain the matrix element, it is very important to obtain h^ t (1; rriQ, tuq') with a procedure 
as similar as possible to the one used to obtain h + (uj;mQ,mQi) for to ^ 1. Thus, we 
get (1; rriQ, itiq/) from the time-component of the ratio of Eq. (B^) with p' = q = 0. 
For degenerate transitions, there is another zero-recoil channel, which corresponds to the 
forward scattering of a meson with one unit of lattice momentum. We do not use the h^ t (1) 
from this channel to determine Zy because it is statistically much noisier than the one 
at zero momentum, and because it does not correspond to a zero-recoil transition in the 
non-degenerate case. 
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Now, to subtract the discretization errors that Zy incorporates, we simply define the 
continuum form factors to be 

m Q , m QI ) ee (1 + /?+(!; m Q , m Q ,)) h ^ ^ m ^ 



h-(u;m Q ,m Qf ) = (1 + m Q , m Q ,)) m<3 ' ^ . (38) 

This definition yields 

h + (uj;m Q} m Q ,) ~ [1 + /? + (w; m Q , m Q ,) + 7 + (w; m Q , m Q ,) - 7 + (l; m Q , m Q ,) 

+d + (u; m Q , m QI ) - d + (l; m Q , m Q ,)]£(u) , (39) 

up to higher-order discretization errors, radiative and power corrections. It is clear from 
Eq. (|39| ) that part of the discretization errors have been subtracted. The subtraction is only 
complete, however, if d + (u) is a constant. For the form factor h~(oj) it is less clear that 
we are subtracting the relevant discretization errors. Indeed, according to the definition of 
Eq. (|38|) the discretization errors in h~ are (d~(u)— d + (l))+0 (a 2 ). However, the assumption 
behind this subtraction is the same as the one made by Lepage, Mackenzie and Kronfeld 
f28fl in their attempt to remove discretization errors by modifying the normalization factors 
which match fermion fields to their continuum counterparts. 

We wish to emphasize here that our subtraction procedure removes non-perturbatively 
all discretization errors which do not break Euclidean invariance and does so to all orders 
in a. Thus, amongst others, all discretization errors which are removed in mean-field theory 
by the procedure of Kronfeld, Lepage and Mackenzie will be removed non-perturbatively by 
our procedure. 

As Eq. ( |3"9"D indicates, in subtracting discretization errors in h + , we also subtract the 
zero-recoil power corrections, 7 + (l), thereby losing the ability to determine them. This is 
not a serious concern in practice because these ought to be small — they are proportional 
to the square of the inverse heavy-quark mass — and therefore difficult to isolate reliably. 
It does mean, however, that even if we can reduce all of our errors to the percent level, we 
will be unable to obtain the zero-recoil power corrections to the form factor Ha-i relevant for 
B — > D*£u decays if we use an analogous subtraction procedure. This is unfortunate because 
these l/m 2 -corrections are one of the dominant theoretical uncertainties in the extraction 
of the CKM matrix element V c b from experimental studies of these decays (see Section [VI H ). 



For obtaining the Isgur-Wise function, however, the fact that our normalization proce- 
dure subtracts these zero-recoil power corrections, which are non-perturbative and difficult 
to quantify, is an advantage. Our hope is that, once these corrections are subtracted, the 
resulting form factor will have smaller power corrections away from zero recoil. 

There is one additional issue surrounding normalization that we wish to address. As 
indicated in the previous subsection, the discretization errors on our three-point functions 
are typically larger for spatial than for temporal channels (see Table jVT|) . Thus, we ought 



to normalize spatial and temporal channels differently. For degenerate transitions, this is 
possible because there is a zero-recoil three-point function which has a non-zero spatial 
component: C^(t;ap = {n/12, 0, 0), 0)q^q. As mentioned above, however, this three-point 
function does not correspond to a zero-recoil decay when Q ^ Q'. We have no zero-recoil 
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three-point function with a non- vanishing spatial component for non-degenerate transitions 
(momenta are quantized on the lattice). So, in order to treat degenerate and non-degenerate 
transitions in the same way, we will normalize h + and h~ as described in Eq. (|38|). 

It is important to note that because h + is obtained from the temporal component of 
Eq. ( pop alone (see end of Section [11 C| ) and is correctly normalized, it does not suffer from 
the possible discrepancy in normalization between temporal and spatial channels. It is 
h~, obtained from both temporal and spatial components, which in fact will absorb this 
discrepancy. For degenerate transitions, where h~ is in principle zero, the values of h~ that 
we obtain are therefore an indication of how large an error this discrepancy can induce 
in the form factors. For non-degenerate transitions, the values of h~ we obtain, though 
contaminated to some extent by discretization errors, can be used to put bounds on the 
physical h~. 



IV. THE FORM FACTORS h+{uj) AND h'(u) 
A. Results at Fixed Light-Quark Mass 



In Table |VII| - Table fX] we present the measurements of h + (uj), h + (u))/(l + (3 + {uj)) 
and h~(oj) which we obtain for all available combinations of the initial and final heavy- 
quark masses for light antiquarks with K q = 0.14144 (Table 0.14226 (Table |VTTI|) 
and 0.14262 (Table |TX| ). In these tables, the first x 2 /d.o. /.-column corresponds to the fit 
which yields h^ t from the temporal component of the ratio W 1 assuming /ij~ t (u) = 0. The 
second x 2 /d.o. /.-column corresponds to the fit which gives h^ at (lv) from both temporal and 
spatial components when holding /ij£ t fixed to its temporal-component value. The number of 
degrees of freedom (d.o.f.) that we quote in this second column depends on the momentum 
channel because the number of non-vanishing equations for hf^ and varies with initial 
and final meson momenta. 

As evidenced by the low values in the first x 2 / d.o.f. column of all three tables, the 
fits which give h^ t from the temporal component of are very good. The fact that 
the values in the second x 2 / d.o.f. column of these tables are generally larger may be due 
to the fact that spatial and temporal components of our three-point functions may have 



different discretization errors, as mentioned in Section |II Q . When we fit these components 
simultaneously to the asymptotic form of Eq. fl25|) while holding fixed, we are not 
fitting to a form which takes into account these discrepancies and consequently obtain a 
larger x 2 1 d.o.f.. As discussed in Section [11 C| , however, this fitting strategy is the only one 
that guarantees that h + does not suffer significantly from discretization errors. 

Given the number of different mass combinations and momentum channels we have, 
our results for h + (u)/(l + (3 + {uj)) are remarkably consistent. Keeping the light-quark mass 
fixed we find that for recoils u which are approximately the same, the values of h + {uj)/{l + 
(3 + {u)) are equal within errors even when they are obtained from different momentum and/or 
heavy-quark mass combinations. This supports the validity of our procedure and is also an 
indication that the radiative corrections obtained using Neubert's results |7| are accurate. 
The fact that h + {uj)/{\ + (3 + (u)) does not appear to depend strongly on the mass of the 
heavy quarks is also an indication that the coefficients of the corrections proportional to 
inverse powers of the heavy-quark masses are not very large (see Section [V]). 
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There are two momentum combinations on which we wish to comment. The first is 
p = (7r/12a, 0, 0) to p' = (7r/12a, 0,0) which, for degenerate transitions, has zero recoil. 
For such transitions, current conservation requires that h + {{ji /12a, 0, 0) — > (it /12a, 0, 0)) 
equal 1. We find values of h + ((n /12a, 0, 0) — > (7r/12a, 0, 0)) which are just barely consistent 
with 1 at the level of 1 a for K q = 0.14144. The situation deteriorates when the mass 



of the light quark decreases (see Table [VIII] and Table |IX]) . Since for given quark masses 
h + ((7r /12a, 0, 0) — > (it /12a, 0, 0)) is extracted from a single three-point function (the one with 
p' = (7r/12a, 0, 0) and q = (0,0,0)), it is much more susceptible to statistical fluctuations 
than most other values of h + which are obtained from averages of three-point functions 
over many equivalent momentum combinations. To show that this slight discrepancy is 
statistical, we consider two measures of h + {{ji /12a, 0, 0) — > (it /12a, 0, 0)) which use the 
same three-point function and normalization. The first is 

fc+((»/12M,0) - (x/12.,0,0)) = (0,0,0) ^ (0,0,0)) 

Zj, // ( Ko ;f I = 4;( I r/12a,0,0) -> (tt/120,0,0)) 

with Zy defined in Section |111 A] . The second is the expression above multiplied by the 
ratio C-i(tf \ (ir/12a, 0, 0))/ 7r/12a) where 7r/12a) is the average of the six p = 

(±7r/12a, 0, 0), (0, ±7r/12a, 0) and (0, 0, ±7r/12a) two-point functions!. Using the values of 
Zy given in Table |VI], the first procedure gives h + ((n/12a, 0, 0) — > (7r/12a, 0, 0)) equal to 
1 to within 1%, even when the mass of the light spectator antiquark is reduced, while the 
second procedure gives results very much in line with the rather low results of Table [VII] , 
Table [VIII| and Table fX[ The reason why the first procedure is more precise is explained in 



footnote |7]. Moreover, that the results given by the second procedure agree better with our 
standard procedure for obtaining h + should not be too surprising, as the latter also makes 
use of average two-point functions. 

The second small inconsistency we wish to comment on is the one arising from the com- 
parison of h + ((0, 0, 0) -> (vr/12a, 0, 0); m Q , m Q >) with h + ((vr/12a, 0, 0) -> (0, 0, 0); m Q >, m Q ), 
for which u is the same. To check the validity of our results, we have re-analysed our data by 
fitting our three-point functions directly to the asymptotic form given in Eq. (|2T|), fixing the 
energies and wavefunction factors which appear in this asymptotic form to their two-point 
function values, and normalizing the resulting h^iu) according to Eq. fl55p. This procedure 
yields values for h + which are nearly identical to the ones given in Tables [VI I| , [VIII| and [TX| . 



The only values that change significantly compared to the size of their error bars are those 
corresponding to h + ((0,0,0) — > (7r/12a, 0, 0); mq, mgr). In this different way of analyzing 
the data, the values we find for h + ((0,0,0) — > (7r/12a, 0, 0); tuq, raqi) are lower, making 
them nearer the values for h + ((7r/12a, 0, 0) — > (0, 0, 0); rngr, toq). This partial bridging of 
the gap, however, comes at the expense of large x 2 /d.o.f.'s, ranging from 2 to 5. One can 
fix both problems — bridging the gap completely and bringing the ^/d.o.f. down — by 
fitting the time component of our three-point functions to 

C 3 (t;p',qW -> MP^(P' 2 ) e - { E P - Ep , +S E )t -E P , tf ( p V) | Vj (0)|p(p)) (41) 



This ratio of two-point functions should of course be 1 in the limit of infinite statistics. 



17 



with an extra parameter SE, instead of to the form given in Eq. (|2l|). The parameter SE is 
designed to absorb slight statistical differences in the time behavior of two- and three-point 
functions. One would worry about the consistency of adding this extra parameter if it were to 
be large compared to the values of the various energies which enter the exponential factor in 
Eq. ([II]) since it is inconsistent to allow for changes in the energies while holding wavefunction 
factors fixed — the two quantities are extremely correlated — and it is inconsistent to claim 
that Ep — Epi is different for two- and three-point functions but that Epi is the same. 
However, we find values of SE which are on the order of 10 -3 and consistent with zero. 

In addition to reconciling the values for h + ((0, 0, 0) — ► (n/12a, 0, 0); mg, mg') and 
h + ((7r/12a, 0, 0) — > (0, 0, 0); rriQi, mg), this method increases the statistical errors on all val- 
ues of h + (uj) because of the additional freedom introduced by the new parameter. We do 
not use this new fitting method as our main one because of the potential inconsistencies 
mentioned above and because the introduction of the extra parameter SE is difficult to 
generalize sensibly to situations where one simultaneously fits more than one four-vector 
component of a three-point function. 

The results given by all of these different methods of analyzing the data are consis- 
tent within statistical errors. This gives us faith that the results for h + in Table [VIH , 



Table |V111| and Table [FX] are valid representations of our data. The most likely rea- 



son, then, for the slight discrepancy between h + ((0, 0, 0) — ► (n/12a, 0, 0); mg, rriQ>) and 
h + ((n/12a, 0, 0) — > (0, 0, 0); mg<, itlq) is that it arises from the same statistical fluctua- 
tion that yields the low value for h + ((jr/12a, 0, 0) — > (jr/12a, 0, 0); mg, mg). Like the 
three-point function which gives h + ((n/12a, 0, 0) — > (n/12a, 0, 0); mg, mg), the one from 
which h + ((7r/12a, 0, 0) — > (0, 0, 0); mgi, mg) is obtained is not averaged with equivalent 
three-point functions. h + ((0, 0, 0) — > (n/12a, 0, 0); mg, mgi), on the other hand, is ob- 
tained from the average of the six three-point functions corresponding to the transitions 
(0,0,0) -> (±7r/12a,0,0), (0, ±vr/12a, 0), (0, 0, ±vr/12a). 

As mentioned earlier, current conservation requires that h~(u) = for degenerate tran- 
sitions. In order to determine whether our results are consistent with this requirement, 
we must know how large the discretization errors on h~(u) might be. As suggested by 
the results for Z v * (see Table [VTD , there may be discretization errors of the order of 10% 
which cause the spatial components of our three-point functions to be low compared to the 
temporal components. One can easily convince oneself, by considering the set of equations 
corresponding to different components of the vector current in Eq. fl25|) , that such discretiza- 
tion errors would cause |/i~(c<j)| to take on values up to about 0.1. This is indeed what we 
find. Thus, to the level of accuracy with which we can determine \h~(uj)\, we can conclude 
that h~(ou) is consistent with zero for degenerate transitions. 

For non-degenerate transitions, the results we obtain for h~(u) resemble very much those 
found in degenerate transitions. They are consistent with zero at the level of 2<ji. Thus, as 
far as we can resolve, h~(u) is small for all uj, most probably less than about 0.1 to 0.2. 

Using this information, we can now put a bound on the size of the error that we are 



9 h (to) is large with large errors when p = p' = (n/12a, 0, 0), because its coefficient in the equation 
which determines it is v 1 — v = (ir/12)(l/amg — 1/amgi), a small number when mg ~ mg/. 
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making on h + (u) by neglecting the contribution of h~(co) to the temporal component of 
Eq. (p5|). Using the fact that the ratio of velocity factors, r = |t>° — v°'\/(v° + t> 0/ ), is at 
most 0.07, and that h + (uo) is always greater than 0.6, we find that the error we make on 
h + {uj) is at most r max h~ (u) max / h + (u) min ~ 1% to 2%. In most situations, if not all, it will 
be smaller than that. Thus, neglecting the contribution of hr{u) in obtaining h + (oj) is a 
very good approximation indeed. 



B. Chiral Extrapolation 

In the previous section, we determined h + (u) for many different combinations of initial 
and final heavy quarks and for three light antiquarks whose masses straddle that of the 
strange quark. In the present section, we describe the extrapolation of our results for 
h + (uj)/(l + (3 + (uj)) to vanishing light-antiquark mass for which K q = K crit =0.14315(2) ||29|| . 
The chirally extrapolated results are relevant for the study of semi-leptonic decays of heavy- 
light mesons whose light antiquark is a u or a d. These results are summarized in Table [X|. 

The extrapolations are covariant and linear in the improved, bare quark mass, am q = 
am q (l — (l/2)am q ), where am q = (l/2n q — l/2n cr - lt ). We fit h + /(l + (3 + ) and u to the 
forms aj l +(am I q ) + j3 h + and a^am^) + (3 U , respectively. Then, h^ Tit /(l + (3 + ) = (3h+ and 
^crit = Pui- The x 2 1 d.o.f. for these extrapolations are given in columns four and six of 
Table [X|. As evidenced by the small values of these x 2 1 d.o.f '.'s, the extrapolations are for the 
most part very smooth. The only extrapolation which has an anomalously large x 2 /d.o.f. is 
the one for h + ((n/12a, 0, 0) — > (n/12a, 0, 0); rriQ, rriQi). As mentioned in Section [IV A| , even 
though current conservation requires that h + ((it /12a, 0,0) — > (7r/12o, 0, 0); rriQ, vtiqi) = 1 
when Q = Q', our results do not quite satisfy this constraint due to a statistical fluctuation. 
As Tables |V11| , |V111| and [TX] further indicate, this constraint is less and less well-satisfied as 
the mass of the light quark is reduced. The correlated extrapolation appears to correct for 
this downward trend in the data, but does so at the expense of a large x 2 /d.o.f . 

We do not extrapolate h~(uo) because this form factor potentially suffers from rather 



large discretization errors as discusssed in Section II C and is therefore not entirely physical. 



C. Interpolation to the Strange Quark 

In the present section, we describe the interpolation of our results for h + (u) /(l + (3 + (u)) 
in the light-antiquark mass to the mass of the strange quark (k s =0.1419(1) |2S[). The 
interpolated results are relevant for the study of semi-leptonic decays of heavy-light mesons 
which contain a strange antiquark. The results are summarized in Table [XI]. They are 
obtained from the same covariant, linear fits as the chirally-extrapolated results of Section 
|IVB| so that the \ 2 j d.o.f . are the same as in Table [X[ The only difference is that the 
interpolated results are hf/(l + j3 + ) = a h +(am I s ) + @ h + and u s = a u (onij) + (3 U , where ml 
is the improved, bare mass of the strange quark. 
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V. DEPENDENCE OF h+(u) ON HEAVY-QUARK MASS 



Having obtained h + (u)) to good accuracy, we can now attempt to determine the depen- 
dence of this quantity on the masses of the initial and final heavy quarks. For the purpose of 
this study, we have calculated h + {uj)/(l + (3 + {oj)) for additional heavy-quark combinations 
when K q = 0.14144. We concentrate on the results which correspond to our heaviest, light 
antiquark (k 9 = 0.14144) because these results have smaller statistical uncertainties and will 
therefore enable us to resolve the dependence of these results on heavy-quark mass more ac- 
curately. We will assume, in the following, that our findings for K q = 0.14144 provide a good 
description of the behavior on heavy-quark mass of our results for smaller light-antiquark 
masses. That this assumption may be justified is confirmed by the mild dependence of h + (oj) 
on light-antiquark mass (see Section |VT]). 

The first indication that the dependence of h + (uj)/(l + j3 + (uj)) on heavy-quark mass 
must be very weak is shown in Fig. [|. In this figure, we plot together the form factors 
h + {uj)/{l + j3 + (u)) for each of our four Q — > Q, degenerate transitions with n q = 0.14144. 
It is natural to begin looking for small heavy-quark mass effects in this data because its 
normalization is free of uncertainties associated with radiative or power corrections (see 
Section pTCl) . 



The four sets of data lie very much on the same curve. To show this more precisely, we 
fit each set individually to the parametrizations £nr{w) and s£nr(u where 

W«) = ——7 (-(2p 2 - (42) 
uj + 1 V 0J + 1J 

is a parametrization for the Isgur-Wise function suggested by M. Neubert and V. Rieckert 
[IJ. In Eq. (f42|), p 2 = — £'(1). We have introduced the supplementary parameter s to 



m 



absorb possible normalization errors. We summarize our findings in Table [XI]| and plot the 
fit curves in Fig. [|. These results clearly show that the four different data sets are entirely 
compatible and suggest that the dependence of h + (u)/(l + f3 + {uj)) on heavy-quark mass is 
expected, therefore, to be quite small over most of the range of experimentally accessible 
recoils 0. 

Before interpreting this observation, let us quantify this heavy-quark-mass dependence 
more precisely. We will do so under the assumption that this small dependence is due to 
power corrections. We have also tested the assumption that it is due to amg-discretization 
errors but find that this assumption is less well satisfied by our data (i.e. it leads to higher 
X 2 /d.o.f.). According to Eq. (||), we have 



10 Since we are only interested in comparing h + (uj)/(l + (3 + {uj)) for different heavy-quark mass, 
any reasonable parametrization will do. 

n The fact that the values of x 2 /d-o.f. are relatively high for all of these fits is explained after 
Eq. pi) in Section [VI. 
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Now, to leading order in the heavy-quark expansion, 



7 + (cj; m Q , m Q >) = g Q (u, a s (m Q ) ,z) e Q + g Q > (u, a s {m Q ,) , z) e Q > 

+0(e 2 Q ,e Q „e Q e Q ,) , (44) 

where 6q(o = A 4 i4 4 /(2mQ(;)) @, and z = mQ/rriQi. The functions qq and qqi correspond 
to matrix elements of dimension-five operators in the HQET Lagrangian evaluated at order 
O (cqm)- These two functions must be equal when Q = Q'. They must also be equal in the 
absence of radiative corrections as HQET cannot distinguish the flavour of a heavy quark at 
order O (e^o) • 111 the presence of radiative corrections, however, the two functions will have 
different values when Q ^ Q'. The amount by which they differ will be partly governed by 
logarithms of the heavy-quark masses, as indicated by the presence of the running coupling 
constant in the functions' arguments. The way in which gg and gQ> depend on z will also be 
different. Nevertheless, since the difference between gq and gQ' is a difference of radiative 
corrections, it is very small. We will neglect this difference in what follows and assume that 

l + {u;m Q ,m Q ,) = g{uj) {e Q + e Q ,) + £> (4, ej,, e Q e Q /) . (45) 

It is worthwhile noting, at this point, that Luke's theorem requires 

0(1) = . (46) 

To evaluate g(u) we need h + (uj;mQ,mQ/) at a fixed uj for different Q or Q'. Because 
momenta on the lattice are quantized this is difficult to achieve. There is one kinematical 
situation, however, where we have enough measurements of h + (uj) at fixed uj for different 
heavy quarks to determine g(uj). When the momentum of one of the mesons vanishes, 
uj becomes independent of that meson's mass. There are four values of uj for which this 
happens, corresponding to |p| = 7r/12a and |p'| = for kq = 0.121, 0.129, and |p| = and 
|p'| = 7r/12a for kqi = 0.125, 0.133. For each of these four points, we have four measurements 
of h + (uj) corresponding to four different values of the mass of the meson which is at rest. 
We pick one of these four measurements and use it to normalize the remaining three. Thus, 
we construct the ratio: 

R+ , s_ 1 L h + (u-m Ql m QI )/{\ + f3 + (u]m Ql m Q >)) \ 
e Q i V h+(u;m Q i,m Q >)/(l + P+(u;m Q i,m Q >)) ) 
= g(uj)(l-x) + 0(e Q ,e Ql ) , (47) 

with x = rriQi/mQ. Here we have assumed that it is the initial meson which has vanishing 
momentum. We then fit the resulting three data points for R + at fixed u to a straight line 
in x. The slope and intercept of this line is g(ui) (see Eq. (0)). We summarize the details 
of these fits in Table |Xlllf In Fig. we show this data with the corresponding fit (solid 



line) for each one of the four values of uj. The data for R + satisfies the parameterization of 



12 Here, A4144 is the energy carried by the light degrees of freedom when n q = 0.14144. We take 
it to be A4144 = ^ (3(i\4 144 - My) + (M 4144 - M*)) + A x = 0.63 GeV. 
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Eq. (f47|) surprisingly well. One should remember that all power corrections are subtracted 
at uj = 1 by our normalization procedure (see Eq. (^) and ensuing discussion). However, 
this is not a problem if one is interested only in 0(eQ,eQr) power corrections to h + since 
these must vanish at zero recoil according (Eq. (|46|)). 

In Fig. |6| we plot g as a function of uj. g{uj) is consistent with zero over the range of recoils 
uj that we can explore (1 < uj < 1.1). Since g{oo) shows no trend over that range and since 
the functions h + (uj) / (1 +[3 + (oj)) plotted in Fig. |3] exhibit no mass dependence over a range of 
recoils from 1 to 1.4, we conclude that g{uj) ought to remain small (less than about 0.2) over 
the full range of experimentally interesting recoils (1 < uj < 1.55). We believe that these 
results indicate that the 1/mQ-corrections to h + (uj) and the remaining amQ-discretization 
errors in h + (uj) are genuinely small because we explore a non- negligible range of heavy-quark 
masses — from about 1 to 2 GeV. It seems quite unlikely that discretization errors or higher 
order power corrections would cancel the leading power corrections over such a range. 

Because g{ui) appears to be less than about 0.2 over the full range of recoils, we predict 
that power corrections to the form factor h + corresponding to physical B — » Div decays 
must be less than about O (e 2 c ) ~ 3% to O (0.2 x (e b + e c ), ef) ~ 5% to 10% over the full 
range of recoils for m b = 4.80 GeV, m c = 1.45 GeV and A = 0.50 GeV Thi s is 

significantly smaller than the O (e c ) ~ 15% {O (e 2 c ) and O (e&) may each contribute an 
additional 5%) that one may naively have expected. It appears, then, that the protection 
Luke's theorem provides at zero recoil extends over the full range of recoils and that for 
the particular combination h + (uj)/(l + (3 + {uj)) the flavour component of the heavy-quark 
symmetry is well satisfied in the charmed sector. This is in stark contrast with our findings 
for the decay constant, fp, of the pseudoscalar D meson ||12|| . In Ref. [|12[] we find that the 
O (e c ) corrections to the heavy-quark limit prediction for this decay constant are of the order 
of 30%. 

These results for g(u) also mean that our results for h + (uj)/(l + f3 + {u))) are, to a good 
approximation, infinite heavy-quark-mass results. Thus, the functions h + (oj)/(l + (3 + {uj)) 
that we measure are effectively Isgur-Wise functions and we can consistently combine data 
corresponding to different initial and final heavy quarks. This is what we do in the following. 

In principle one could also try to quantify power corrections to h~(oj). In the absence 
of radiative corrections, we find from the results of Ref. that these power corrections are 
given by: 

l~(w; m Q: m Q >) = (1 - 277(0;)) (-e Q + e Q >) , (48) 
where, like g(w) defined in Section [V], rj(u) is a subleading, universal form factor 0. Eq. 



indicates that power corrections proportional to eq are equal and opposite to those propor- 
tional to €q'. This prediction is consistent with the mass dependence we observe in our 
results. However, because our normalization procedure is optimized for determining h + (oj) 



13 We have included, in this estimate, potential higher order corrections that may have been 
subtracted by our normalization procedure. 



14 



Luke's theorem does not constrain r](uj) at u = 1 as it did g{ui). 
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and not h~(u), it is not clear to what extent the mass- dependence due to power corrections 
can be resolved from that due to discretization errors and to higher-order power corrections 
coming from our normalization procedure (see Eq. (B3)). 



VI. DEPENDENCE OF h+(u) ON LIGHT-QUARK MASS: ISGUR-WISE 

FUNCTIONS 

We established in the previous section that, for fixed light- ant iquark mass, we can com- 
bine the results for h + (u) / (1+/5 + (c<j)) corresponding to different values of the initial and final 
heavy-quark mass. We further established that the resulting function is an Isgur-Wise func- 
tion: £ Ui d(co>) when the mass of the light antiquark vanishes; £ s (u;) when the light antiquark 
is given the mass of the strange quark. 

We plot £, u ,d(u) and £ s (u;) in Fig. [7] and Fig. ||, respectively. We fit the corresponding 
data to the parametrizations s£,nr(u), t; NR (u). The parameter s is added to absorb possible 
uncertainties in the normalization of these form factors. Because the parametrization £,nr{&) 
is only one of many possible parametrizations, we also fit our results to s£u n (uj), £u n (u), 
£,quad{u) and s£ quad (u) where 

£, m (u) = l- p 2 {u-l) (49) 
is a simple linear parametrization, and 

W) = 1 - f?(u " 1) + - l) 2 (50) 



is a quadratic parametrization. The parameter c in Eq. fl50[) is, of course, the curvature of 
the Isgur-Wise function at u — 1. We tabulate the results of these different fits in Table |XTV| . 
In this table, we also present the results of performing these fits on the data corresponding 



to K q =0.14144, 0.14226 and 0.14262. 



The fact that the values of x 2 / d.o.f. are relatively high for all of these fits should 
not in itself be taken as an indication that the parametrizations of Eq. (f4"2"D, Q4T3D and 
(p0|) are poor representations of the Isgur-Wise functions. These large values of x 2 /d.o.f. 
are due to the discrepancy that we mentioned in Section IV A| between our measurements 



of h + ((0, 0, 0) -> {ir /12a, 0, 0); m Q , m QI ) and of h + ((vr/12a, 0, 0) -> (0, 0, 0); m Q ,, m Q ). Be- 
cause of this discrepancy, no parameterization can fit our data with a good value of x 2 /d.o.f.. 
X 2 /d.o.f '.'s nevertheless seem to favor the use of the extra parameter s but does not seriously 
discriminate between s^r(u), s£u n (u)) and s £, qua d{u) . We have tried fitting our data to 
yet other parametrizations and of all the fitting functions, s£n n (u) yields the lowest values 
for p 2 . The reason for this is that s£u n (u) is the only parametrization which does not have 
positive curvature. Since s£u n (uj) is in that sense an exception, we will not use it as our 
standard fitting function but because it is a valid parametrization for these Isgur-Wise func- 
tions we will make certain that our results have errors which encompass the values it gives 
for the slope. Furthermore, since both s£,nr(u) and s£ qua< i(ui) give nearly identical fits (see 
Figs. and ^|), we will use s£nr(u) as our standard in the following because it has one less 
parameter and yields better x 2 /d.o.f.'s. 



23 



Having already argued in Section [V] that, with the normalization that we have adopted, 
the remaining mass- dependent discretization errors are small, we turn to momentum- 
dependent lattice artefacts. To quantify these momentum-dependent discretization errors 
we resort to the following procedure. We fit the data for £ U) d(u) and £ a (w) for fixed ini- 
tial and final meson momentum and all heavy-quark combinations, to the parametrization 
s£nr(uj). The variation in the results of fits to these different momentum sets should give 
us some indication of how large these momentum-dependent lattice artefacts are. Some of 
this variation, of course, may be due to statistical fluctuations of the sort we mentioned in 
Section [V~&[ 

We summarize the results of the fits to the different momentum sets in Table IX VI It is 



reassuring that the value of s for the case (0, 0, 0) — > (n/12a } 0, 0) is very close to 1, because 
the corresponding data are our best points. They are the points for which our normalization 
procedure is optimal because they are obtained from three-point functions which are much 
more correlated with the three-point function which yields the normalization factor h^ t (1)0. 
Furthermore, these data have the smallest statistical errors and should have the smallest 
discretization errors, because the momenta of the incoming and outgoing mesons are less 
than or equal to the initial and final momenta of other momentum sets. 

To accommodate the spread in the values in the slope parameters p\ d and p 2 s correspond- 
ing to £ u ,d(w) and £ s (w), we assign systematic errors to these parameters which encompass 
all the central values given in Table |XV| . The central value and statistical errors that we 
quote are given by fitting s£nr(u) to all momentum sets put together (see Table |XdV| ). 
Thus, our final results for the slope at u = 1 are 

PM = 0.9^(stat.)^(syst.) (51) 

for i u4 {u) and 

^ = 1.2^(stat.)l?(syst.) (52) 

for £ g (w). Even though the exact values of these slope parameters are slightly different if 
different parametrizations for the Isgur-Wise functions are used, these differences are well 
within our error bars. 



In Table XVI we compare our predictions for the slope of the Isgur-Wise functions £ u a 



and £ s with those of other authors. We find that our predictions for p 2 lie safely above 



the lower bound of Bjorken and below the upper bound of de Rafael and Taron [32 
Our results for p 2 also agree with the lattice result of Bernard et al. || obtained with 
Wilson fermions for a light spectator antiquark with mass m q ~ m s , although the details 
and systematics of the two calculations are different. The authors of Ref. |§ do not quote 
a value of p 2 u d for vanishing light-quark mass. 



5 For readers familiar with the methods used to calculate three-point functions, the reason why 
the three-point functions corresponding to (0,0,0) — > (ir/12a, 0, 0) and (0,0,0) — > (0,0,0) are 
strongly correlated is because they are built up from the same exponentiated propagator. Indeed, 
the initial momentum in our notation is the momentum of the exponentiated propagator. 



24 



Also for comparison, we quote an average experimental value for the slope of the Isgur- 

and 



Wise function compiled by Neubert [33] from very recent results of the ALEPH [fl 
CLEO [|T7| Collaborations as well as older data from the ARGUS Collaboration [31 



PM (e ^.) = 0.87(12)(20) 



(53) 



where the second error is theoretical and accounts for the uncertainty in the size of l/m c 
corrections [Q. Eqs. (|51| ) and ( |53"D agree remarkably well. 



As can be inferred from Eqs. (|5lf) and fl52|) and from Table XIV| our results are compatible 
with the statement that p 2 is constant with light-quark mass, possibly decreasing slightly 
as this mass decreases. Such a decrease is consistent with one's intuition that it is more 
difficult to make the light degrees of freedom recoil, the heavier these degees of freedom are. 



Furthermore, H_0gaasen and Sadzikowski |35| find a decrease in slope which is very similar 
to the trend we observe in the central value of p 2 when we include the extra parameter s in 
our fits. In fact, our predictions for p 2 itself are in excellent agreement with theirs. Their 
prediction is based on an improved bag model calculation and is an extension of earlier work 
by Sadzikowski and Zalewski ||36|| . A similar decrease in slope with spectator quark mass 

are 



and p s 



is observed by Close and Wambach J57J] though the values they quote for p 2 d 
slightly larger than the ones we find. 

To test the robustness of our predictions for p 2 , we have explored many different pro- 
cedures for obtaining h + (u), two of which we have already described in Section IV A . To 
obtain h + (u) for degenerate transitions, we have in addition tried normalizing our lattice re- 
sults for ht &t {u) by ((7r/12a,0,0) (n/Ua, 0, 0)) instead of by h^ t , ((0,0,0) (0,0,0)) 
(see Eq. (|38|)). When fitted to the s£(u;) parametrizations, the results obtained using all of 
these methods give very similar values for the slope parameter p 2 . They only differ in the 
value of s they predict, i.e. in their overall normalization. Thus, we are quite confident that 
our predictions for the slope are reliable but believe that it is important to allow for the 
extra normalization parameter s. 



VII. EXTRACTION OF V CB 

In Section [TV] we obtained h + (oj) for a variety of P — > P' transitions where P (P') is 
a pseudoscalar meson composed of a heavy quark Q {Q') and a light antiquark q. In our 
study, both Q and Q' are quarks with masses around that of the charm quark. In Section |V], 
however, we showed that our results for h + {u)/{l + (3 + {uj)) are independent of heavy-quark 
mass for masses around the charm quark mass or larger. This means, modulo the issue of 
power corrections at zero-recoil, that our results can be used to describe not only P — > P' 
transitions with Q{Q') ~ c but B q — > D q decays as well, where the subscript q labels the 
flavour of the light antiquark. In Section[VT], we studied the dependence of h + {u) / {1+ j3 + {u)) 
on the mass of the light, spectator quark, m q , and obtained results for ^^(cu) and ^(cu). All 
of this means that our result for once multiplied by (1 + /3^ c ), is in fact the form factor 
h + relevant for B u ^ —>■ D u ^ transitions, while (1 + /3^ c )^ s is the form factor h + relevant for 
B s — > D s transitions. 
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Now, the differential decay width for B — > D£z/ is, in the limit of zero lepton mass, 
dT(B (s) ^D {s) e,) _ W (Mfi(s) + Md j, M , y _ 1)3/2 



(54) 



V / Af S(s) +Af D(a) V / 



So, in principle, we could obtain \V c b\ by comparing our theoretical prediction for 
dr ( B ( s ')~^ D ( s ) eu ^ t an experimental measurement of this rate. A major problem with this 

approach, however, is that the rate '^^M^Elil^L j s helicity suppressed, as evidenced by the 
factor (to 2 — 1) 3//2 , so that it is very difficult to get accurate experimental measurements close 
to oo — 1 where the predictions of HQET are most reliable. Another problem with obtaining 
\Vcb\ from JB( S ) — > D^tv decays is that one must know hr(uo) to better accuracy than is 
given by our calculation: an error of 0.1 on h~ leads to an uncertainty of about 10% in the 
rate. We should mention, however, that Neubert || has estimated h~ using perturbation 
theory and sum rules in HQET and has found that its magnitude does not exceed 0.04 over 
the whole range of recoils oo. If this is true, its contribution to the rate of Eq. ([54]) should 
not exceed 4%. 

We have not exhausted the predictions of heavy-quark symmetry. We have yet to exploit 
the spin component of this symmetry. Using a combination of the spin and flavour symmetry, 
we can relate our predictions for ~ h + (oo)/(l + f3 + (uo)) to the form factors required to 
describe B^ — > Df^iv decays. These form factors are defined by || 

(D{ s) (p>,e)\crb\B {s) (p)) # , 

— — — = = in (oo)e^ p e,.v'vg , 

M B( Mp* K J v a P 

W (s) 

(u + l)e^h M (u)-e* •v(y' t h Aa (w)+v'>'h A *(u])) , (55) 



(D* {s) (p' 7 e)\crr^b\B (s) (p)) 



M B{S) M D 



where v = p/Mb (s) and v' = p'/Mp* . In the heavy-quark limit, these four form factors 
can be expressed in terms of the single Isgur-Wise function, £, u> d( s )- There are, of course, 
radiative and power corrections to these heavy-quark symmetry predictions. Thus, one has 

h\u) = (a 1 + (3\u) + Y M)£ M(S) M , (56) 

with i = V, Ai, A 2 , A 3 and 

a* = a M = a v = 1 , 

a M = . (57) 



Luke's theorem |l(J further guarantees that, at zero recoil, h Al is free of O (e&, e c ) corrections, 
i.e 7" 4l (l) ~ 0(e 2 ) e 2 c ). Because h Al is the only form factor to contribute to the differential 
decay rate for B^ — ► DJ.lv decays at zero recoil, Luke's theorem implies that the leading 
non-perturbative corrections to this rate must be small at u — 1. More precisely, in the 
limit of zero lepton mass, 
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^Ml u {M B(a) -M Du f [l + ^ (1 



M 2 B(s) -2uM B(s) M DU +M^ 



(M Bw - M D 



K{u) 



(58) 



where P Al (l) = —0.01 0] and K(u>) incorporates the radiative corrections, /^(u;), away 
from to = 1, the non-perturbative power corrections, 7 j4i (cj), and the contributions of the 
three other form factors to the rate. From what we have said above, it should be clear that 
K(l) = l + O (ef, e 2 c ). One can also show [[| that in the limit of exact heavy-quark symmetry, 
K(uj) = 1 for all lu. Moreover, since we have factored out £ u ,d(s)(w) in the expression for 
the rate, the contributions of the three other form factors will be normalized by C«,d(s)- This 
means that K(u) is a collection of radiative and power corrections (see Eqs. (|56|) and fl57|) ), 
many of which are kinematically suppressed: deviations of K(ui) from 1 ought to remain 
small. Neubert estimates [f|, using perturbation theory and sum rules in HQET, that K(u) 
may reduce the slope parameter, p 2 by 0.09 which corresponds to an enhancement of the 
rate by about 10% at maximum recoil and by a smaller amount for smaller uj. However, we 
cannot estimate yet how the physical K(u) deviates from its value in the heavy-quark limit 
from our lattice calculation. For that we need to study Qq(0~) — > Q'q(l~) decays, which we 
are currently analyzing. We also need to determine the 1/m^-corrections to h Al (l), which 
as discussed in Section |111 C| , we cannot get with a procedure analogous to the one presented 
in this paper. Hence, we will assume that K(to) = K(l) for all uj which is a reasonable 
assumption given the size of our errors on the slope parameter p 2 . We can then use our 
lattice determination of £ Ut d{w) to extract V c b from the experimentally measured differential 
decay rate for B — ► D*iv (B s —>■ D*£u has not yet been measured). This analysis differs 
from Neubert 's extraction of \ V c b\ |Tj| in that we fix the cu-dependence of the differential decay 
rate using our calculation instead of fitting it from experiment. This enables us not only to 
extract \V c b\ with no free parameters, but also to check the validity of non-perturbative QCD 
against experiment. We find that the cu-dependence predicted by our calculation agrees very 
well with the results of the ALEPH |T§ and CLEO |T7| collaborations. 

In Fig. ^ Fig. [K] and Fig. |ll| we show least-x 2 -fits to experimental data for |V^,|(1 + 
f3 M {l))K(uj)i u4 (uj) from ALEPH [|T§, ARGUS g§ and CLEO |TJ, respectively. The only 
parameter is \V c b\. The slope of the Isgur-Wise function is constrained to the value given 
by our lattice calculation (see Eq. ([H])) and the functional form for the Isgur-Wise function 
that is used is of Eq. (E3) E3. The results of these fits are summarized in Table PC VII . 



Our results favor ALEPH and CLEO data over that of ARGUS. Using the data from CLEO, 
for instance, we find 



w, 



cb\ 



0.037 



-2 + 4 
2-1 



0.99 



1 + 



1 + h/r. 



(59) 



16 As can be seen from Fig. ^ other parametrizations give very similar curves when fit to our results 
for Therefore, results for \V c b\ obtained with these other parametrizations will be well within 
our quoted error bars. 
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where 8i/ m s are the power corrections proportional to in K(l) which have been the 

subject of much controversy of late []33l ,[38ll . 

For comparison, wepresent recent experimental predictions for |V^,|i^(l) obtained from 
a linear fit to the dataliZI 



V cb \ K(l) 



'1+^(1)' 
0.99 



0.0351 ±0.0019 ±0.0020 
0.0385 ±0.0044 ±0.0035 
0.0392 ±0.0043 ±0.0025 



CLEO g7| 
ALEPH 
ARGUS 



16] 



(60) 



34 1 



where the first error is statistical and the second systematic. These results have been rescaled 



by Neubert using the new lifetime values t b q = 1.61(8) ps and t b + = 1.65(7) ps |j40 
These new lifetimes reduce | V c b\ (l+(3 Al (l)) K{1) by approximately 1%. Our results compare 
very well with these experimental measurements, especially in the case of the CLEO result. 
This is due to the fact that our Isgur-Wise function has an u dependence which agrees very 
well with that of the CLEO data. 



VIII. EXCLUSIVE DECAY RATES 



Having determined the Isgur-Wise functions £ Uj d and £ s , 
ing ratios. For — > D^iu decays, all we have to do is integrate Eq. ( 
the results by the B( s ) meson lifetime. We approximate h + in Eq. ([54]) by 
p 2 given by Eq. (|5T]) or Eq. (52) depending on whether the light antiquark is a u, d or an 



we can evaluate exclusive branch- 
[54] ) and multiply 
{1 + /3 + )Znr with 



s (see Section |V11| ). We neglect the O (eb,e c ) power corrections to h + since they appear to 
be small (see Section |V|) but add a 10% error to account for possible higher order power 
corrections. We further neglect the contribution of h~ in accordance with Neubert's findings 
that this form factor is smaller than 0.04 over the whole range of u (see discussion after 
Eq. d54D ) but add 20% to our errors since an \h~(u)\ ~ 0.2£(u;) is consistent with radiative 
corrections and order of magnitude estimates of power corrections. 



The branching ratios for B, 



D 



(a)* 



decays are equally simple to obtain. Here it 



D* 



{Ml {s) +Ml u )/2M B{sy ^ i 

decays. As discussed 
= 1 which leads to an 



D {a) £u 



is Eq. (pq) that we must integrate over the range 1 < u < 

The Isgur-Wise functions used are the same as for B^ a \ — > 
after Eq. (]58|), we assume K{uS) = K(l). We further assume K{uS) - 
uncertainty of the order of O (2e^) ~ 5 — 10% in the branching ratios 

We summarize our results for B (B( s ) ~^ D^iv J and B (B( s ) 
The first set of errors is obtained by adding our lattice statistical and systematic errors in 
quadrature. The second set of errors corresponds to the uncertainty due to deviations from 
the heavy-quark limit. For comparison, we list the experimentally measured values for these 
branching ratios. Agreement with our predictions is very satisfactory. 



in Table XVIII 



D 



Finally, we give a prediction for the ratio of the rates for B^ — > Df^iv and for B 



(s) 



In this ratio, the factors of \V c b\ cancel and lifetimes do not appear. This ratio is 



thus a purely theoretical prediction. We find 
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The ARGUS result has been corrected for the new D branching fractions [39|. 
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T(B -> Dlv 
and 



I(Z? D '" ] =3.2^(lat.)±1.0(hqs) (61) 



T(B s ^D s £u] 



3.3IJ(lat.) ± l.O(hqs) . (62) 



where the first set of errors was obtained by adding our lattice statistical and systematic 
errors in quadrature and the second set of errors, denoted by "hqs" , quantifies the uncertainty 
due to neglected power and radiative corrections. For comparison, the experimental result 
for T Y^§^De.D) is 2.1(1) [[|lj. Though low compared to our prediction, this result is consistent 
with ours within errors. 



IX. CONCLUSIONS 

We have presented an extensive study of semi-leptonic B^ — > Dt s )£v decays where we 
evaluate the matrix element, (D\c^f^b\B) , for many different values of mb and m c around the 
physical charm mass and three values of the light antiquark mass around that of the strange. 
Because the charm quark has a bare mass which is almost 1/3 the inverse lattice spacing, 
mass-dependent discretization errors are a problem that we must contend with. To reduce 
these errors we use an O (a)-improved quark action in which the leading such errors are 
no longer O (amg) but rather O (a s arriQ, (arriQ) 2 ). This reduces discretization errors from 
O (40%) to O (5 — 15%) at the mass of the charm. To reduce them even further we describe, 
in Section |111 C| , a procedure for subtracting them at least partially. Only those discretization 
errors which have the same dependence on uj as h + (oj) will be fully subtracted. We believe, 
however, that the observation in Section [V| of /i + (u;)'s lack of dependence on heavy-quark 
mass indicates that a fairly large proportion of discretization errors are eliminated with our 
procedure. 

The fact that we obtain h + (uj) and hr{ui) for many values of the initial and final heavy- 
quark masses enables us to study their heavy-quark-mass dependence. We find that the 
residual dependence of h + /(l + f3 + {uj)) on the heavy-quark mass is consistent with zero. 
Given our errors, we conclude that power corrections to the form factor h + for physical 
B — > D transitions are less than 10%. This is much smaller than the 25% corrections one 
is entitled to expect for form factors not protected by Luke's theorem. It is also in stark 
contrast with our findings for the decay constant, fp, of the pseudoscalar D meson 



In Ref. [12] we find that the O (e c ) corrections to the heavy-quark limit prediction for this 
constant are on the order of 30%. Thus, it appears that the protection from O (Aqcd/^c) 
effects that Luke's theorem provides at zero recoil extends to some extent over the full 
range of experimentally accessible uj. Our results for h + {uj)/ (1 + (3 + {uj)) are then, to a good 
approximation, the corresponding Isgur-Wise function. 

Having obtained the Isgur-Wise function from h + (u)/(l + (3 + {oj)) for three values of the 
mass of the light, spectator antiquark, we can study its dependence on light-quark mass. 
Interpolating the light antiquark to the strange, we obtain an Isgur-Wise function relevant 
for B s — > D^W decays which has a slope — £' s = 1.2^2(stat.) ^(syst.) at zero recoil when 



fit to a parametrization proposed by Neubert and Rieckert pOf . Extrapolating to a massless 
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light antiquark yields an Isgur-Wise function relevant for B — > D^*'W decays. This function 
has a slope —£' ud = O.Ql^stat.) * 2 (syst.) at zero recoil. We observe a slight decrease in the 
magnitude of the central value of the slope as the mass of the light antiquark is reduced 
in accordance with one's understanding that more massive degrees of freedom have more 
inertia. Given the errors, however, the significance of this observation is limited. 

We also use these functions, in conjunction with heavy-quark effective theory, to extract 
V c b from the experimentally measured B — > D*£u decay rate. Our procedure for extracting 
\Vcb\ differs from that proposed by Neubert |l5j in that we fix the cu-dependence of the 
differential decay rate using our calculation instead of fitting it from experiment. This 
enables us not only to extract \V c b\ with no free parameters, but also to check the validity 
of non-perturbative QCD against experiment. We find that the a;-dependence predicted 
by our calculation agrees very well with the results of the ALEPH |L6| and CLEO |Ti 
collaborations. Using the data from CLEO, for instance, we find 

1^1 = 0.0371}^^ 

where 5i/ m § is the power correction proportional to l/m 2 c at zero recoil and (3 Al (l), the 
relevant radiative correction. Here, the first set of errors is due to experimental uncertainties, 
the second due to statistical errors and the third to systematic errors in our evaluation of the 
Isgur-Wise function. We also use our Isgur-Wise functions and heavy-quark effective theory 
to calculate branching ratios for B^ — > D^iv and £>( s ) — > D1 8 \tv decays. Agreement with 
experiment is very good. Finally, we compute the following ratios of rates: 

|^l = 3.2^(lat.)±1.0(hqs) (63) 

and 

T(B S -> D*£u 




T(B s ^D s £u] 



3.3ti(lat.) ±1.0(hqs) , (64) 



where the first set of errors was obtained by adding our lattice statistical and systematic 
errors in quadrature and the second set of errors, denoted by "hqs" , quantifies the uncertainty 
due to neglected power and radiative corrections. In these ratios, the factors of \V c b\ cancel 
and .B-meson lifetimes are absent: they are purely theoretical predictions. 

We are currently extending our study to the matrix elements relevant for B^ — > DJ s \li> 
decays. This will enable us not only to check our predictions for the various Isgur-Wise 
functions but also to test the heavy-quark spin symmetry. We are also undertaking a study 
of semi-leptonic A& — ► A c and — > S c decays, where the Ab( c ) is a J p = l/2 + baryon 
composed of a 6(c) quark and two light quarks coupled to spin and isospin and the H 6 ( c \, 
another J p = l/2 + baryon composed of a b(c) quark and two light quarks but this time 
with spin 0, isospin 1/2 and strangeness —1. That study should provide many interesting 
phenomenological predictions which are at the limit of current experimental knowledge as 
well as many tests of the heavy-quark symmetry. Finally, we are planning to repeat these 
studies on lattices with different lattice spacing in order to remove discretization errors in a 
more systematic way. 



30 



ACKNOWLEDGMENTS 

L. L. would like to thank Southampton University's theory group for its hospitality and 
financial support during much of this project. We would also like to thank G. Martinelli 
and members of the APE Collabortaion for stimulating discussions. This research was sup- 
ported by the UK Science and Engineering Research Council under grants GR/G 32779 and 
GR/J 21347, by the European Union under contract CHRX-CT92-0051, by the University 
of Edinburgh and by Meiko Limited. We are grateful to Edinburgh University Computing 
Service and, in particular, to Mike Brown for maintaining service on the Meiko i860 Com- 
puting Surface. CTS (Senior Fellow), DGR (Advanced Fellow) and NMH acknowledge the 
support of the Particle Physics and Astronomy Research Council. JN acknowledges the Eu- 
ropean Union for for their support through the award of a Postdoctoral Fellowship, contract 
No. CHBICT920066. 



31 



REFERENCES 



[1] For a comprehensive review of B decays see, for example, "B Decays", ed. S. Stone, 2 nd 

Edition (World Scientific, Singapore, 1994). 
[2] S. Nussinov and W. Wetzel, Phys. Rev. D36, 130 (1987); M. B. Voloshin and M. A. Shif- 

man, Sov. J. Nucl. Phys. 45, 292 (1987); 47, 511 (1988). 
[3] N. Isgur and M. B. Wise, Phys. Lett. B232, 113 (1989); Phys. Lett. B237, 527 (1990). 
[4] H. Georgi, Phys. Lett. B240, 447 (1990). 

[5] UKQCD Collaboration, S. P. Booth et al., Phys. Rev. Lett. 72, 462 (1994). 
[6] CW. Bernard, Y. Shen, A. Soni, Phys. Lett. B317, 164 (1993). 
[7] M. Neubert, Phys. Rev. D46, 2212 (1992). 

[8] A. F. Falk, H. Georgi, B. Grinstein and M. B. Wise, Nucl. Phys. B343, 1 (1990). 
[9] M. Neubert, Phys. Rept. 245, 259 (1994). 
[10] M.E. Luke, Phys. Lett. B252, 447 (1990). 

[11] UKQCD Collaboration, C.R. Allton et al, Nucl. Phys. B407, 331 (1993). 
[12] UKQCD Collaboration, R. M. Baxter et al., Phys. Rev. D49, 1594 (1994). 
[13] B. Sheikholeslami and R. Wohlert, Nucl. Phys. B259, 572 (1985). 

[14] G. Heatlie, G. Martinelli, C. Pittori, G. C. Rossi and C. T. Sachrajda, Nucl. Phys. B352, 
266 (1991). 

[15] M. Neubert, Phys. Lett. B264, 455 (1991). 

[16] I. Scott (ALEPH Collaboration), to appear in the proceedings of the 27th International 
Conference on High Energy Physics, Glasgow, Scotland, July 1994. 

[17] T. Browder (CLEO Collaboration), to appear in the proceedings of the 27th International 
Conference on High Energy Physics, Glasgow, Scotland, July 1994; B. Barish et al. (CLEO 
Collaboration), Phys. Rev. D51, 1014 (1995). 

[18] UKQCD Collaboration, C.R. Allton et al., Phys. Rev. D47, 5128 (1993). 

[19] CW. Bernard and A. Soni in Quantum Fields on the Computer, ed. M. Creutz, (World 
Scientific, Singapore, 1992). A.S. Kronfeld and P.B. Mackenzie, "Progress in QCD using 
Lattice Gauge Theory", Annual Review of Nuclear and Particle Science 43, 793 (1993). 

[20] G. Martinelli, Nucl. Phys. B (Proc. Suppl.) 9, 139 (1989). 

[21] B.Efron, SIAM Review 21, 460 (1979). 

[22] M. Beneke and V.M. Braun, Nucl. Phys. B426, 301 (1994). 

[23] I. Bigi et al, Phys. Lett. B339, 160 (1994). 

[24] UKQCD Collaboration, D.S. Henty, R.D. Kenway, B.J. Pendleton and J.I. Skullerud, 
Edinburgh Preprint 94/545, Phys. Rev. D (in press). 

[25] A.Borrelli, R. Frezotti, E. Gabrielli and C. Pittori, Nucl. Phys. B409, 382 (1993). 

[26] G.P Lepage and P.B. Mackenzie, Phys. Rev. D48, 2250 (1993). 

[27] D. Henty (UKQCD Collaboration), Nucl. Phys. B (Proc. Suppl.) 34, 468 (1994). 

[28] G. P. Lepage, Nucl. Phys. B (Proc. Suppl.) 26, 45 (1992). G. P. Lepage and P. B. Macken- 
zie, Phys. Rev. D48, 2250 (1993); A. S. Kronfeld, Nucl. Phys. B (Proc. Suppl.) 30, 
445 (1993); P. Mackenzie, |hep-lat/92120Tl talk given at DPF92, Batavia, Illinois, 10-14 
November 1992; Nucl. Phys. B (Proc. Suppl.) 30, 35 (1993). 

[29] UKQCD Collaboration, C.R. Allton et al, Phys. Rev. D49, 474 (1994). 

[30] M. Neubert and V. Rieckert, Nucl. Phys. B382, 97 (1992). 

[31] J. D. Bjorken, SLAC report SLAC-PUB-5278 (1990) and "La Thuile Rencontres", 583 
(1990). 



32 



E. de Rafael and J. Taron, Phys. Rev. D50, 373 (1994). 
M. Neubert, Phys. Lett. B338, 84 (1994). 

H. Albrecht et al. (ARGUS Collaboration), Z. Phys. C 57, 533 (1993). 
H. H0gaasen and M. Sadzikowski, Z. Phys. C64, 427 (1994). 
M. Sadzikowski and K. Zalewski, Z. Phys. C59, 677 (1993). 

F. E. Close and A. Wambach, Rutherford Appleton Laboratory preprint RAL-94-041, 
April 1994, |hep-ph/94053T4] . 



M. Shifman, TPI-MINN-94/31-T, July 1994, |hep-ph/9409358 



R. Patterson, to appear in the proceedings of the 27th International Conference on High 
Energy Physics, Glasgow, Scotland, July 1994. 

F. DeJongh, to appear in the proceedings of the 27th International Conference on High 
Energy Physics, Glasgow, Scotland, July 1994. 

"Semi-leptonic B Decays", S. Stone, in "B Decays", ed. S. Stone, 2 nd Edition (World 

Scientific, Singapore, 1994). 

S. Narison, Phys. Lett. B325, 197 (1994). 

M. B. Voloshin, Phys. Rev. D46, 3062 (1992). 

B. Blok and M. Shifman, Phys. Rev. D47, 2949 (1993). 

J. L. Rosner, Phys. Rev. D42, 3732 (1990). 

G. Burdman, Phys. Lett. B284, 133 (1992). 

Y. B. Dai, C. S. Huang and H. Y. Jin, Z. Phys. C56, 707 (1992). 

W. Venus, private communication. 

Particle Data Group, Phys. Rev. D45, SI (1992). 



33 



FIGURES 

FIG. 1. The ratio R ' (i) 3 up to constant factors, vs. t for the case where the initial meson 
has momentum (0,0,0) and the final meson, momentum (7r/12a, 0, 0). Here, the initial and final 
heavy-quark hopping parameters are kq = kqi =0.129 while the light-quark hopping parameter is 



K q =0.14144. The solid line is obtained from our fit of R°(t) to the asymptotic form of Eq. (25). 



The dashed lines indicate the errors of this fit. 

FIG. 2. Values of Zy^ as functions of mga. The solid lines represent fits to quadratic functions 
of mqa for the data at the two different values of (5. We have also plotted the light-quark values 
of Zyf given in Eq. ( |30|) but have not included them in the fit. 

FIG. 3. Values of Zy ] (kq) obtained from the simulation at j3 = 6.2 at different momenta and 
Lorentz indices. The three curves are quadratic fits to the three sets of data. 

FIG. 4. h + {uj)/(l + (3 + {uj)) vs. u for all four elastic scattering reactions: kq = kqi = 0.121, 
0.129, 0.133. The light-quark hopping parameter is fixed to K q = 0.14144. The curves are obtained 
by fitting each heavy-quark, data set to s^nr(^)- The data points as well as the curves correspond- 
ing to different heavy quarks are really indistinguisible. The 1/mg corrections to h + (u) / (1+ f3 + (u)) 
cannot therefore be very large. (See text for details.) 

FIG. 5. R + vs. x = mQ'/rriQ at fixed oj for four values of u. The solid lines are obtained by 
fitting these results to the parameterization given in Eq. (|47| ) and the dotted lines represent errors. 
The slope and intercept of this line is the subleading form factor g(u) (Eqs. ( |45| ) and (p3|)). The 
light-quark hopping parameter is K q = 0.14144. 

FIG. 6. The subleading form factor g{oS) (Eqs. fl45| ) and (f43"D). The light-quark hopping pa- 
rameter is K q = 0.14144. 

FIG. 7. £ U) d(u;) = ^ + ( Ll -')/(l +/9 + ( Ct; )) vs - u f° r K q = K crit- The different symbols correspond to 
different values of initial and final heavy-quark mass. The solid curve is obtained by fitting these 
results to s£at_r(w) while the dashed curve corresponds to a fit to s£n n (uj) and the dotted curve 
corresponds to a fit to s£ qua d(u>). The value of p 2 shown on the plot is the one given in Eq. (|5l|). 



FIG. 8. £ s (u;) = h + {uo)/{\ + (3 + (uj)) vs. u for n q = k s . The different symbols correspond to 
different values of initial and final heavy-quark mass. The solid curve is obtained by fitting these 
results to s^nr(uj) while the dashed curve corresponds to a fit to s£n n (uj) and the dotted curve 
corresponds to a fit to s^ qua( i(uj). The value of p 2 shown on the plot is the one given in Eq. (|52"1). 
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FIG. 9. Least-x 2 -fit to experimental data for |V cfe |(l + j3 Al (1)) K{u) £(u) from ALEPH @ 



assuming K(tu) = K(l). In this fit, the slope of the Isgur-Wise function is constrained to the value 



given by our lattice calculation (see Eq. (51)) and the functional form for the Isgur-Wise function 
that is used is (,nr of Eq. ([42|). The first set of errors on | V c b\ is due to experimental uncertainties, 
the second set of errors results from the lattice statistical errors on p 2 , and the third, from the 
lattice systematic errors on p 2 . The experimental points were obtained from a measurement of the 
rate dB{B — > D*lv) / duj. Also shown are our appropriately scaled, chirally-extrapolated results 
(octagons). 



FIG. 10. Same fit as in Fig. |9| but for experimental data from the ARGUS Collaboration [34|. 



FIG. 11. Same fit as in Fig. |9| but for experimental data from the CLEO Collaboration [17|. 
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TABLES 



TABLE I. Physical heavy-quark masses corresponding to different values the heavy-quark hop- 
ping parameter, kq. They are obtained from the corresponding chirally-extrapolated pseudoscalar 
and vector meson masses, as described in Eq. (ffif). For completeness, we also tabulate the chi- 
rally-extrapolated meson masses in lattice units (a" 1 ss 2.7 GeV @). These masses were obtained 
by covariant linear extrapolation of the masses Mp and My obtained at three values of the light 
antiquark hopping parameter: K q = 0.14144, 0.14226, 0.14262. The pseudscalar meson masses 



were computed as described in Section II D, while the vector meson masses were obtained as in 
PJ, with a fitting range 11 < t < 23. 





Ml 


My 


m Q (GeV) 


0.121 


0.874 +3 


0.896 +4 


1.90 


0.125 


0.773 +3 


0.799 +3 


1.64 


0.129 


0.665 +3 


0.696 +4 


1.36 



0.133 0.547 ±3 0.588 +5 1.06 



TABLE II. Wavefunction factors, Z 2 , and energies, Ep for our heavy-light, pseudoscalar 
mesons and for two values of momentum, |p|. The energies are quoted in lattice units 
(a -1 ~ 2.7 GeV [Q). The x 2 /d-o.f. for the fits which give these results are all on the order 



of 1. 





K Q 




0.121 






0.125 






0.129 






0.133 






IpI 


Z 2 




Ep 




z 2 




E P 




z 2 




Ep 








Ep 




0.14144 




vr/12a 


17.9 
12.3 


+ 6 
-5 
+ 5 
-5 


0.924 
0.958 


+ 2 
-2 

+ 3 
-2 


16.3 
11.4 


f 5 
-5 
f 4 
-5 


0.823 
0.861 


+ 2 
-2 

+ 3 
-2 


14.5 
10.3 


f 4 
-4 
f 4 
-5 


0.716 
0.760 


+ 2 

-2 

+ 3 

-2 


12.4 ; 
9.0 t 


f 4 
-4 
-4 
4 


0.600 
0.653 


+ 2 
-2 

+ 3 
-3 


0.14226 




vr/12a 


15.5 
10.5 


+ 5 
-5 
+ 5 
-5 


0.901 
0.937 


+ 3 
-2 
+ 3 
-4 


14.2 
9.7 j 


f 5 
-5 
-4 
-4 


0.800 
0.840 


+ 3 
-2 
+ 3 
-3 


12.7 

8.8 ^ 


f 4 
-4 
-4 
-4 


0.692 
0.739 


+ 3 
-2 
+ 3 
-3 


10.8 ; 
7.7 1 


f 5 
-3 
-4 
3 


0.575 
0.631 


+ 3 
-2 
+ 4 
-3 


0.14262 





14.7 


+ 8 
-7 


0.892 


+ 4 

-4 


13.5 


f 7 
-6 


0.791 


+ 3 
-3 


12.0 


f 6 
-5 


0.683 


+ 3 
-3 


10.3 


f 4 
-4 


0.565 


+ 3 
-2 




vr/12a 


9.8 j 


-7 

-5 


0.928 


+ 5 
-4 


9.1 1 


-6 
-4 


0.832 


+ 4 
-4 


8.3 j 


-5 
-4 


0.730 


+ 5 
-4 


7.3 1 


-5 
4 


0.623 


+ 4 
-3 



36 



TABLE III. (3 + (lo) vs. lo for all combinations of initial and final heavy-quark mass. 



KQ — ► KQI 


1.0 


1.1 


1.2 


L.3 


1.4 


0.121 -► 0.121 





-0.025 


-0.047 


-0.068 


-0.088 


0.121 -v 0.125 


0.017 


-0.008 


-0.030 


-0.051 


-0.071 


0.121 -► 0.129 


0.037 


0.013 


-0.009 


-0.030 


-0.050 


0.121 -► 0.133 


0.063 


0.040 


0.018 


-0.002 


-0.022 


0.125 -► 0.125 





-0.023 


-0.045 


-0.065 


-0.085 


0.125 -► 0.129 


0.024 


0.001 


-0.021 


-0.041 


-0.060 


0.125 -► 0.133 


0.055 


0.033 


0.012 


-0.008 


-0.027 


0.129 -► 0.129 





-0.022 


-0.042 


-0.061 


-0.079 


0.129 -► 0.133 


0.039 


0.017 


-0.003 


-0.022 


-0.039 


0.133 -► 0.133 





-0.019 


-0.038 


-0.055 


-0.071 


TABLE IV. vs. w 


for all combinations of initial and final heavy-quark mass. 








1.0 


1.1 


1.2 


1.3 


1.4 


0.121 -► 0.121 

















0.121 -► 0.125 


0.001 


0.000 


-0.001 


-0.001 


-0.002 


0.121 -» 0.129 


-0.003 


-0.004 


-0.005 


-0.006 


-0.007 


0.121 -► 0.133 


-0.014 


-0.016 


-0.017 


-0.019 


-0.021 


0.125 0.125 

















0.125 -> 0.129 


0.000 


-0.001 


-0.001 


-0.002 


-0.003 


0.125 0.133 


-0.008 


-0.009 


-0.011 


-0.012 


-0.014 


0.129 0.129 

















0.129 -»• 0.133 


-0.001 


-0.002 


-0.004 


-0.005 


-0.006 


0.133 -»■ 0.133 
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TABLE V. Values of the effective normalisation constant as a function of the improved 

bare mass of the heavy quark. The value of K q is 0.14144 at (3 = 6.2 and 0.144 at (3 = 6.0. 





/3 = 6.2 








(3 = 6.0 






w} Q a 




) 




m} Q a 




0.133 


0.231 


0.8913 t 


2 
1 


0.129 


0.344 


0.920 ± \ 


0.129 


0.310 


0.9177 + 


3 
2 


0.125 


0.405 


0.945 + \ 


0.125 


0.379 


0.9428 t 


4 
2 


0.120 


0.464 


0.973 + 2 


0.121 


0.435 


0.9659 t 


5 
3 









TABLE VI. Values of Zy* for different choices of the Lorentz index /j,, momenta p, and light 



quark masses (given by 


K g ) from the simulation at /3 = 6.2. 






H and p 






Zy 3 








K q = 0.14144 


K q = 0.14226 


^ = 0.14262 


H = 4, p = 


0.133 


0.8913 t I 






fi = 4, p = 


0.129 


0.9177 t 2 


0.9168 t 4 


0.9165 t e 


/i = 4, p = 


0.125 


0.9428 + 2 






H = 4, p = 


0.121 


0.9659 + 3 


0.9656 + jj 


0.9658 + ^ 


/i = 4,p= (tt/12,0,0) 


0.133 


0.8976 1 X ° 






^ = 4,p= (tt/12,0,0) 


0.129 


0.9248 t ? 


0.9242 + ^ 


0.9240 1 24 


^ = 4,p= (tt/12,0,0) 


0.125 


0.9498 1 I 






^ = 4,p= (tt/12,0,0) 


0.121 


0.9729 1 J 


0.9734 t H 


0.9746 t 25 


A* = 1,2?= (tt/12,0,0) 
H = 1,P= (tt/12,0,0) 
H = l,p= (tt/12,0,0) 
H = l,P= (tt/12,0,0) 


0.133 
0.129 
0.125 
0.121 


0.949 + % 
0.994 + S 
1.042 t 67 
1.084 + ™ 


0.982 + H 
1.089 + g 


0.924 ^ 
1.059 3SS 
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TABLE VII. Results for h + {uj), h + (uo)/(l + (3 + {uj)) and h~{uj) obtained with the fitting pro- 
cedure described in Section |II C| . The light-quark hopping parameter is fixed to K q = 0.14144 and 
all heavy-quark mass combinations are presented. Only transitions with initial and final meson 
momenta less or equal to 7r/(12a) are included. 



1 } 




UJ 




h + (uj) 




/l+( 




v 2 Id o f 


hr 




\ 2 Id n f 










= 0.121 






= 121 k„ = 


0.14144 








(0,0,0) 


(1,0,0) 


1.037 


+ 1 

— 1 


0.95 t 


1 
1 




0.96 + } 

— 1 


3.6/2 


0.12 


+ 2 
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TABLE VIII. Results for h + {uj), h + (uj)/(l + (i + {uj)) and hr(u) obtained with the fitting pro- 
cedure described in Section |II C| . The light-quark hopping parameter is fixed to K q = 0.14226 and 
all heavy-quark mass combinations are presented. Only transitions with initial and final meson 
momenta less or equal to 7r/(12a) are included. 
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TABLE IX. Results for h + (uj), h + (u)/(l +P + (u)) and h~(u) obtained with the fitting proce- 
dure described in Section |II C| . The light-quark hopping parameter is fixed to K q = 0.14262 and 
all heavy-quark mass combinations are presented. Only transitions with initial and final meson 
momenta less or equal to n/(12a) are included. 
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TABLE X. Results for to and h + (to)/(l + (3 + (lo)), for n q = K crit =0.14315(2), obtained from 
covariant, linear extrapolations of the results for n q =0.14144, 0.14226, 0.14262. All heavy-quark 
mass combinations are presented. The first y 2 /d.o.f. column corresponds to the uj extrapolation; 
the second, to the extrapolation of h + (u))/(l + [3 + {lo)). 
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TABLE XI. Results for u and h + (u)/(l + /3 + (w)), for K q = k s =0.1419(1), obtained from 
covariant, linear interpolations of the results for K q =0.14144, 0.14226, 0.14262. All heavy-quark 
mass combinations are presented. The x 2 /d.o.f. are the same as for the chiral extrapolations (see 
Table g). 
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(1,0,0) 




1.007 tl 


0.95 + ° 


(1,0,0) 


(0,0,0) 




1.039 +2 


0.91 t\ 


(1,0,0) 


(0,1,0) 




1.136 + _t 


0.81 


(1,0,0) 


(-1,0,0) 




1.264 +g 


0.75 +3 




Kq = 0.129 - 


-+ K Q' 


— . 133 5 K/q — K/g 




(0,0,0) 


(1,0,0) 




1.094 +3 


0.90 +} 


(1,0,0) 


(1,0,0) 




0.998 +g 


n ^ . -4- 

0.94 +5 


(1,0,0) 


(0,0,0) 




1.065 +2 


0.87 I? 


(1,0,0) 


(0,1,0) 




1.165 +g 


0.78 +2 


(1,0,0) 


(-1,0,0) 




1.332 +7 


0.70 +2 



TABLE XII. Results of fits of our data for h + (u) /(l + [3 + (u>)) to the parametrizations s £jv_r(u;) 
and £,nr{u) described in the text. The four kq = k'q(= 0.121,0.125,0.129,0.133) transitions with 
K q = 0.14144 are considered in turn. Only transitions with initial and final meson momenta less 
or equal to (n/12a) are included. 















£,nr{u) 


KQ = k'q 








X 2 /d.o.f. 


P 2 


X 2 /d.o.f. 


0.121 


1.4 


+ 3 
-3 


0.99 +} 


12.5/3 


1.6 + 2 3 


13.0/4 


0.125 


1.4 


+ 2 
-3 


0.99 +} 


13.6/3 


1.6 + _l 


14.1/4 


0.129 


1.4 


+ 1 
-2 


0.99 t\ 


13.5/3 


1.5 + _l 


13.9/4 


0.133 


1.4 


+ 1 
-2 


0.99 t\ 


11.1/3 


1.4 


11.5/4 
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TABLE XIII. Power corrections to h + {u) for four values of uj when n q = 0.14144. See text for 
definition of R + and g{uj) 



pt 






PQ 


PQ' U 


X 




1 


0.121 


0.121 
0.125 
0.129 

U. loo 


(1,0,0) 


(0,0,0) 1.037 + } 


1.00 
1.16 
1.40 

1 7Q 




-0.015 til 
-0.031 +23 

-U.UD/ _ 35 








n(, A — 

g^) — 


07^ a/ 2 /Hr»f — 1 /0 
U.U(0 _g]^ Willi ^ /QUI — U.l/Z 






2 


0.121 
0.125 
0.129 

U. loo 


0.125 


(0,0,0) 


(1,0,0) 1.047 t\ 


1.00 
1.16 
1.40 

1 7Q 




0.002 +7 
0.021 +^ 
n n^i 

U.UOl 








n(, A — 


0/11 ^~ 2 ^ Titri + ln a/ 2 /^^f fi A /9 

— U.U41 _29 Willi \ /Qui — U.4/Z 






3 


0.129 


0.121 
0.125 
0.129 

U. loo 


(1,0,0) 


(0,0,0) 1.062 +2 


1.00 
1.16 
1.40 

1 7Q 




-0.006 t 13 9 
-0.049 +25 
-u.uoz 44 








q(uj) = 


0.083 +gg with x 2 /dof=0.8/2 






4 


0.121 
0.125 
0.129 
0.133 


0.133 


(0,0,0) 


(1,0,0) 1.088 +2 


1.00 
1.16 
1.40 
1.79 




0.010 t\l 
0.013 tH 
0.003 tf 4 








$M = 


-0.025 with x 2 /dof=0.5/2 
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TABLE XIV. Results of fits of our data for h + (oj) / (1 + j3 + (lo)) to the parametrizations £at_r(u;), 
^Hni^) an d ^quadi 1 ^) with and without the additional parameter s, as described in the text. For 
fixed K q , all heavy-quark mass combinations are used. Only transitions with initial and final meson 
momenta less or equal to 7r/(12a) are included. Here K cr i t =0.14315(2) and k s =0. 1419(1). 







P 


2 




S 


X 2 /d.o.f. 




P 2 


X 2 /d.o.f. 


0.14144 




1.3 


+ 2 

-1 




0.98 


109/38 




1.5 ^ 


121/39 


K s 




1.2 


+ 2 
-2 




0.98 t\ 


95/38 




1.4 ^ 


106/39 


0.14226 




1.0 


+ 2 
-3 




0.96 +i 


113/38 




1.4 ±1 


140/39 


0.14262 




0.7 


+ 3 
-3 




0.94 +2 


100/38 




1.4 ^ 


134/39 


^crit 




0.9 


+ 2 
-3 




0.96 ±2 


69/38 




1.3 % 


88/39 














S6in(w) 










Kq 




P 


2 




S 


X 2 /d.o.f. 




P 2 


X 2 /d.o.f. 


0.14144 




1.0 


+ 1 




0.97 ±\ 


111/38 




1.3 +1 


159/39 


K s 




0.9 


+ 1 




0.97 


97/38 




1.2 


139/39 


0.14226 




0.8 


+ 1 
-2 




0.96 t\ 


114/38 




1.2 ^ 


170/39 


0.14262 




0.6 


+ 2 
-2 




0.93 +2 


100/38 




1-1 IS 


155/39 


^crit 




0.7 


+ 1 
_ 2 




0.95 I2 


71/38 




1.1 ^ 


111/39 












Stquadi^) 










Kq 


P 


2 




C 


s 


X 2 /d.o.f. 


P 2 


C 


X 2 /d.o.f. 


0.14144 


1.2 


+ 2 
-2 


1.6 


+ 1.2 
-1.3 


0.98 li 


108/37 


1.6 ^ 


3.9 t\l 


115/38 


Kg 


1.2 


+ 2 
-2 


2.0 


+ 1.2 
-1.5 


0.98 +2 


94/37 


1.5 ^ 


4.0 ^ 


99/38 


0.14226 


1.0 


+ 2 

-3 


1.4 


+ 1.3 
-1.7 


0.96 t\ 


113/37 


1.6 + _l 


5-1 t\i 


125/38 


0.14262 


0.7 


+ 4 
-4 


1.1 


+ 2.2 
-2.6 


0.94 +^ 


100/37 


1.7 U 


6.9 


114/38 


^crit 


1.0 


+ 3 
-3 


2.3 


+ 1.5 
-2.0 


0.97 +3 


69/37 


1.5 I? 


5.0 t£5 


74/38 
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TABLE XV. The different momentum sets for £ u ,d(w) and £ s (u;) are fit to the parametrization 
s£nr(u>) (Eq. (|4"2"D), A momentum set comprises all combinations of initial and final heavy quarks 
with fixed initial and final momenta. These different fits are used to estimate remaining systematics 
(see text). 



















6( 


u) 




p 


p' 


P 2 




s 


X 2 /d.o.f. 


P 2 




S 




x7d.o./. 


(0,0,0) 


(1,0,0) 


1.1 1 


-3 
-2 


1.01 


0.3/6 


1.2 1 


2 
1 


1.00 


+ 2 
- 1 


0.1/6 


(1,0,0) 


(0,0,0) 


1.3 1 


-5 
- 1 


0.93 + J 


0.8/6 


1.4 1 


3 
3 


0.95 


+ 3 
-2 


0.7/6 


(1,0,0) 


(0,1,0) 


1.2 j 


-4 
-3 


0.95 tjj 


0.1/6 


1.4 1 


2 
2 


0.96 


+ 4 

-3 


0.1/6 


(1,0,0) 


(-1,0,0) 


0.7 j 


-4 

-2 


0.95 ^ 


0.2/6 


1.1 1 


2 
2 


0.95 


+ 7 

™5 


0.4/6 



TABLE XVI. Comparion of our lattice results for — £' ud (l) and —£' s (l) to the theoretical pre- 
dictions of various authors. 



Reference 




-O 1 ) 




UKQCD 




0.9ig(stat.)lJ(^yst.) 


1.2 + 2 (stat.)+ 2 (syst.) 


Bernard, Shen and Soni || 




1.24(26) (stat.) (26) (syst.) 


de Rafael and Taron ]32j 


p 2 < 1.42 




Close and Wambach [jsj]] 


1.40 


1.64 


Narison [|| 


] 


1.00(2) 




Neubert § 




0.66(5) 




Voloshin ||| 


1.4(3) 




Bjorken |31 


i p 2 > 0.25 


Blok and Shifman Q 


0.35 < p 2 < 1.15 




H0gaasen and Sadzikowski [35] 


0.98 


1.135 


Rosner p5| 




1.59(43) 




Burdman jl 


1 


1.08(10) 




Dai, Huang and Jin 47j 


1.05(20) 





TABLE XVII. Results for \ V cb \ from a fit of \V cb \(l + (3 Ai (1))K(uj)£ nr (lj) to experimental data 
with £,nr(u) fixed by our lattice computation (i.e. p 2 is given by Eq. (|5l|)) and K(lj) = 1. The 
experimental data are obtained from the differential branching ratio for B — > D*tv decays. In the 
| V^t, | column, the first set of errors is due to experimental uncertainties, the second set of errors 
results from the lattice statistical errors on p 2 , and the third, from the lattice systematic errors on 



Experiment 


M 1+ f 9 9 (1) )( 1+ <w 


X 2 /d.o.f. 


ALEPH 
ARGUS 
CLEO II 


0.042(2) + 2 t\ 
0.033(2) +2 -i 
0.037(1) + 2 t{ 


3.0/5 
9.9/7 
4.5/6 
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TABLE XVIII. Our predictions for various branching ratios compared to the experimentally 
measured values for these ratios. Our results are obtained assuming \V c b\ =0.038 |Q, rgo =1.53ps, 



t S o =1.54ps @, M so = 5.28 GeV, M B o = 5.38 GeV, M B+ = 1.87 GeV, M B + = 1.97 GeV, 
M B * + = 2.01 GeV and Mg*+ =2.11 GeV p9 |. Our errors are explained in the text. We only 



consider here semi-leptonic B° and B® decays because the experimental data for charged B meson 
decays are much less precise. The quoted experimental numbers were taken from Ref. [41]. 





B -► Dlv 


B 8 -» D s lu 


B -» D*lD 


B s -> D*Ju 


UKQCD 


1.5 +4±0.3 


1.3 +2± - 3 


4.8 +9±0.5 


4.4 ^5±0.4 


ARGUS 


2.1±0.7±0.6 




4.7±0.6±0.6 




CLEO I 


1.8±0.6±0.3 




4.1±0.5±0.7 




CLEO II 






4.50±0.44±0.44 
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"i 1 1 1 1 1 1 1 r 

0| 



1.5 



<D|V U |D> 



I I i ^* 



1 i 



1.0 



0.5 



ic q =/c Q ,=0. 129, fc q =Q. 14144 



0.0 



J I I L 



J I I L 







10 



20 



1.4 i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



1.2 



1.0 



0.8 




o Zf (0,0,0) M=0 
x Zf (1,0,0) M=1 
o Zf (1,0,0) M=0 



3 i i i i i i i i i i i i i i i i i i i i i 

0.0 0.2 0.4 0.6 0. 

amj, 




1.0 



1.1 

CO 



1.2 



i 1 1 r 



i 1 1 r 



i 1 1 r 



0.06 



^3 



< 
+ 



0.04 



o UKQCD 
ALEPH 



NIX 
✓ IN 




0.02 



o 

> 



„2 n q+2+4 
P =0.9_3_ 2 



|V cb | = 0.042333 



0.00 



J I I L 



1.0 



J I I L 



J I I L 



1.2 



1.4 



1.6 



CO 



i 1 1 r 



i 1 1 r 



i 1 1 r 



0.06 



^3 



o UKQCD 
ARGUS 



✓ IN 



< 
+ 



o 
> 



0.04 



0.02 




^2 n q+2+4 
P =0.9_3_ 2 



|V cb | = 0.033^? 



0.00 



J I I L 



1.0 



J I I L 



J I I L 



1.2 



1.4 



1.6 



CO 



i 1 1 r 



i 1 1 r 



i 1 1 r 



0.06 



^3 



o UKQCD 
CLEO II 



✓ IN 



< 
+ 



o 
> 



0.04 




0.02 



^2 n q+2+4 
P =0.9_3_ 2 



|V cb | = 0.037^ 



0.00 



J I I L 



1.0 



J I I L 



J I I L 



1.2 



1.4 



1.6 



CO 



